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ABSTRACT 

Discussion  of  the  fundamental  elements  and  theory  related  to 
atmospheric  transmittance  Is  presented.  A llne-by-llne  transmission 
computer  program  Is  developed  that  utilizes  the  combined  Doppler- 
Lorentz  (Voigt)  line  broadening  function.  In  addition  a rapid 

algorithm  to  evaluate  the  Voigt  function  with  a maximum  relative 

/ 0,  00  O 

A 

error  of  about  one  part  In  10  Is  described  and  a software  package 
that  processes  the  absorption  line  parameters  necessary  to  calculate 
transmittance  Is  given.  The  results  of  transmittance  calculations 
for  seven  channels  In  the  fifteen  micrometer  band,  corresponding 
to  the  seven  High  Resolution  Infrared  Radiation  Sounder  (HIRS) 
channels,  are  furnished  and  the  procedure  to  compute  a band  averaged 
transmission  discussed.  It  Is  concluded  that  the  Voigt  algorithm 
developed  here  Is  an  excellent  computational  procedure  and  the 
resulting  transmittance  program  correctly  calculates  atmospheric 
transmission. 
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CHAPTER  I 
INTRODUCTION 

The  thin  film  of  gas  that  clings  to  the  surface  of  the  earth 
performs  many  functions,  and  has  a correspondingly  large  number  of 
Important  properties.  To  a meteorologist  the  atmosphere  is  basi- 
cally a sink  for  solar  radiation,  and  its  most  interesting  charac- 
teristic is  the  ability  to  transform  that  energy  into  wind  systems 
and  rainfall  averages.  This  phenomenon  of  radiative  transfer 
through  the  atmosphere  is  too  complex,  in  general,  for  a closed 
solution.  However,  if  the  scope  of  investigation  is  narrowed, 
mathematical  models  can  be  developed  that  are  useful.  One  such 
application,  that  of  energy  transfer  theory,  is  the  retrieval  of 

the  vertical  temperature  structure  of  the  atmosphere  from  remote 

23 

satellite  measurements  [Kaplan  ] . This  requires  not  only  accurate 
measurements  and  a stable  mathematical  inversion  scheme,  but  also 
very  accurate  transmittance  functions.  This  transmittance  is  easily 
calculated  if  atmospheric  absorption  is  known.  The  most  accurate 
method  to  predict  atmospheric  absorption  is  to  add  the  contribution 
of  each  molecular  species,  line-by-line.  These  line-by-line 
calculations  are  very  time  consuming  even  with  today's  high  speed 
computers.  This  document  is  concerned  with  molecular  absorption  and 
the  mathematical  functions  that  may  reduce  the  time  required  for  a 
given  line  calculation.  It  is  organized  into  seven  chapters  as 
outlined  below. 

The  equations  of  atmospheric  transmittance  are  reviewed  in 


Chapter  two.  Additionally  the  composition  of  the  atmosphere  and 
some  of  Its  optical  properties  are  examined.  The  constituents  of 
interest  and  the  processes  by  which  they  absorb  energy  are  touched 
upon.  The  concept  of  an  absorbing  line  is  Introduced  and  the 
phenomenom  of  line  broadening  Is  discussed. 

Chapter  three  defines  atmospheric  absorption  and  transmittance. 
Then  the  transmittance  differential  equation  this  work  Is  to  solve 
Is  developed.  A relation  for  the  combined  effects  of  Lorentz  and 
Doppler  spectral  line  broadening  Is  obtained,  and  finally  the 
frequency  bands  of  Interest  for  temperature  retrieval  from 
satellite  soundings  are  defined  and  the  molecular  absorbers  within 
the  band  identified. 

In  Chapter  four  begins  the  major  contribution  of  this  work. 
Three  methods  of  evaluating  the  Voigt  function  are  presented.  The 
results  of  tests  for  accuracy  and  speed  conducted  for  each  of  the 
methods  are  tabulated  and  discussed. 

In  Chapter  five  the  data  used  in  transmittance  calculations 
and  the  units  of  all  equations  are  covered.  Some  remarks  on  the 
derivation  of  absorption  line  parameters  are  given  and  the  computer 
software  used  to  handle  the  data  is  developed.  The  units  of  the 
data  and  necessary  conversion  factors  are  also  presented. 

Chapter  six  ties  together  the  work  of  the  other  chapters  and 
presents  the  end  product  of  this  research,  the  transmittance  program. 

Finally  there  is  discussion  of  the  results  of  this  work  and 
some  conclusions  are  drawn.  Where  appropriate  the  work  of  a chapter 
is  supported  in  an  appendix. 
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CHAPTER  II 


FUNDAMENTALS  OF  ATMOSPHERIC  ABSORPTION 


2.1  Preface 

This  chapter  discusses  the  composition  of  the  atmosphere  and 
molecular  energy  absorption.  The  atmosphere  Is  presented  first, 
followed  by  an  outline  of  molecular  energy  transitions.  Then  spectral 
line  broadening  Is  treated  and  finally  the  Important  points  of  the 
chapter  are  summarized. 

2.2  The  Atmosphere 

The  gases  that  compose  our  atmosphere  are  listed  In  Table  2.2-1 
18 

[Hudson  ] . Those  who  relative  proportions  are  nearly  constant  up  to 
altitudes  In  excess  of  80  kilometers  are  called  the  permanent  constit- 
uents. The  atmosphere  contains  several  other  gases,  called  the 
variable  constituents;  their  amounts  varying  with  temperature,  alti- 
tude and  location.  Chief  among  them  is  water  vapor,  which  may  consti- 
tute as  much  as  2 percent  of  a humid  atmosphere  at  sea  level.  The 
amount  of  water  vapor  decreases  rapidly  with  altitude.  Another 
variable  constituent  ozone,  is  seldom  observed  at  sea  level.  The  amount 
of  ozone  In  the  atmosphere  Increases  with  altitude  to  a maximum  at 
about  20  kilometers  and  then  decreases  at  higher  altitudes.  Fig.  2.2-1 

depicts  the  variance  of  water  vapor  with  altitude  and  season,  while 

32 

Fig.  2.2-2  reflects  the  correlation  for  ozone  [McClatchey  ]. 


Constituent 


Formula 


Percent  by  Volume 


Nitrogen 

«2 

78.1 

Oxygen 

°2 

20.9 

Argon 

Ar 

0.934 

Carbon  Dioxide 

CO2 

0.032 

--3 

Neon 

Ne 

1.81-10  ^ 

Hel lum 

He 

5.24-10"^ 

Methane 

CH, 

4 

2.0-10 

Krypton 

Kr 

1.14-10  ^ 

Nitrous  Oxide 

N2O 

5.0-10  ^ 

Hydrogen 

«2 

5.0-10  ^ 

Carbon  Monoxide 

CO 

2.0-10  ^ 

Xenon 

Xe 

9.0-10  ^ 

Ozone 

O3 

Varies 

Water  Vapor 

H^O 

Varies 

18 

TABLE  2.2-1.  Composition  of  the  Atmosphere  [Hudson  ]. 
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1.  TROPICAL 

2.  MIDLATITUDE  SUMMER 

3.  MIDLATITUDE  WINTER 


0 4 8 12  16  20 

WATER  VAPOR  (gm/m*) 
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Figure  2.2-1.  Variation  of  Water  Vapor  [McClatchey  ]. 
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In  addition  to  the  gases  listed,  the  atmosphere  may  contain  a 
variety  of  particles  that  may  be  of  concern.  Depending  on  the  nature 
of  the  Investigation  It  could  be  necessary  to  consider  salt  from 
ocean  spray,  fine  dust  from  the  earth's  surface,  or  carbon  particles 
resulting  from  combustion.  Near  cities  or  Industrial  complexes 
atmospheric  pollutants  such  as  ammonia,  hydrogen  sulfide,  and  sulfur 
dioxide  are  present  In  significant  amounts. 

2.2.1  Atmospheric  Temperature 

Having  considered  the  composition  of  the  atmosphere  we  now 
discuss  Its  thermal  properties.  The  earth's  atmosphere  Is  divided 
Into  four  regions  primarily  on  the  basis  of  temperature  gradients, 
as  shown  In  Fig.  2.2-3.  Typically  the  temperature  decreases  with 
altitude  (termed  a positive  lapse  rate)  In  the  troposphere  passing 
through  a minimum  at  the  tropopause,  and  then  rises  with  Increasing 
altitude  In  the  stratosphere.  A region  In  which  temperature  Increases 
with  altitude  Is  termed  an  inversion,  or  Inversion  layer.  At  the 
stratopause  the  temperature  passes  through  a more  or  less  well  defined 
maximum,  then  decreases  through  the  mesosphere  to  a second  minimum  at 
the  mesopause.  Above  the  mesopause.  In  the  thermosphere,  the  kinetic 
temperature  again  Increases  with  altitude,  reaching  a maximum  value 
at  a height  of  several  thousand  kilometers.  Above  the  thermosphere, 
and  demarcated  from  It  by  the  thermopause.  Is  essentially  an  Isothermal 
region,  the  exosphere.  At  any  particular  time  or  place  the  true  tem- 
perature structure  may  differ  appreciably  from  Fig.  2.2-3,  but  in 


ALTITUDE  (km) 


TEMPERATURE  (‘K) 


Figure  2.2-3.  Atmospheric  Nomenclature 
(Adapted  from  [45]) 


PRESSURE  (mb) 


general,  the  curve  shown  is  at  least  qualitatively  correct 
[McEwan  and  Phillips^^] . 

Various  regions  or  layers  of  the  atmosphere  have  properties 
which  are  sufficiently  characteristic  to  justify  their  being  given 
special  names,  most  notable  the  ionosphere.  The  ionosphere  itself 
is  divided  into  three  regions  according  to  altitude  and  electron 
density. 

2.2.2  Atmospheric  Pressure 

The  variation  of  atmospheric  pressure  with  altitude  is  governed 

by  the  hydrostatic  equation,  which  is  derived  with  the  aid  of  Fig. 

34 

2.2-4  [McEwan  and  Phillips  ].  The  pressure  is  P at  an  altitude  Z, 
and  P - AP  at  altitude  Z + AZ. 


Figure  2.2-4.  Variation  of  Pressure  with  Altitude. 
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The  difference  In  pressure  is  due  to  the  weight  of  the  gas  in  the 
column  of  unit  cross-sectional  area  and  length  AZ.  Hence  in  the 
limit 

§;  “ -Pg.  (2.2-1) 

where  g is  the  acceleration  due  to  gravity  and  p is  the  gas  density. 
For  air  we  use  the  equation  of  state  for  an  ideal  gas, 

P = nkT.  (2.2-2) 

Dividing  Eq.  2.2-1  by  Eq.  2.2-2  yields 


^ = _£8_dz  = -^ 

P nlcT  H 


(2.2-3) 
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where  H * , is  termed  the  scale  height  [Whitten  and  Poppoff  ) . 

pg 

Asstnnlng  the  scale  height  to  be  constant,  Eq.  2.2-3  can  be  Integrated 
to  obtain 


(2.2-4) 


which  shows  the  expected  exponential  decrease  of  pressure  with  height 
in  a region  of  constant  g and  T.  The  left  scale  of  Fig.  2.2—3  depicts 
the  decrease  of  atmospheric  pressure  with  altitude. 


2.3  Molecular  Energy  Transitions 


In  first  approximation,  the  energy  (E)  of  an  isolated  molecule 
can  be  presented  in  the  form 


E..  + E + 

trans  vib 


E + E , 
ro  el 


(2.3-1) 


where  is  the  translational  motion,  which  depends  on  velocity 

and  may  assume  any  value,  E^^^  is  the  energy  of  the  electrons,  and 
E ..  and  E are  the  vibrational  and  rotational  energies.  The  last 
three  terms  in  Eq.  2.3-1  are  quantized  [Weldner  and  Sells^^]  and 
take  discrete  values  onxy,  the  values  being  specified  by  one  or  more 
quantum  numbers.  Any  combination  of  quantum  numbers  defines  an 
energy  state  or  quantum  state,  or  term.  Radiation  is  absorbed  or 
emitted  when  a transition  takes  place  from  one  energy  state  to  another. 
The  frequency  (v)  of  the  absorbed  or  emitted  quantum  is  given  by 
Planck's  relation, 

AE  = hv,  (2.3-2) 


where  h is  Planck's  constant.  Equation  2.3-1  is  valid  only  when  the 
interaction  of  different  molecular  motions  is  disregarded.  In  general, 
the  energy  of  a molecule  is 


trans 


E ^.  + E + E , + 
vib  ro  el 


E , - + E , + E 

el-vib  el-ro  vlb-ro 


(2.3-3) 
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where  the  last  three  terms  represent  the  interaction  [Zuev  ].  Thus, 
the  most  general  transition  Involves  simultaneous  changes  of  elec- 
tronic, vibrational,  and  rotational  energy. 
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Tfie  electronic,  vibrational,  and  rotational  energies  of  a mole- 
cule differ  in  order  of  magnitude.  The  energy  of  an  electronic 
transition  is  of  the  order  of  several  electron  volts,  a vibrational 
energy  transition  is  reckoned  In  the  tenths  or  hundredths  of  an  elec- 
tron volt,  and  a rotational  energy  transition  is  in  the  thousandths 
or  ten-thousandths  of  an  electron  volt.  But,  since  a molecule  has 
electronic,  vibrational,  and  rotational  energies,  and  these  energies 
change  simultaneously  on  emission  or  absorption,  the  electronic  and 
vibrational  spectra  are  not  evidenced  in  pure  form.  Depending  on 
the  frequency  range  of  Interest  we  deal  either  with  an  electronic- 
vibratlonal-rotatlonal  molecular  spectrum,  or  with  a vibrational- 
rotational  spectrum,  or  with  purely  a rotational  spectrum.  These 
spectra  are  referred  to  as  the  electronic,  vibrational,  and  rota- 
tional spectra  [Goody^^],  respectively. 

Each  electron  state  has  a corresponding  grid  of  vibrational 
energy  levels  characterized  by  certain  values  of  the  vibrational 
quantum  number  v.  Each  vibrational  level  has  a corresponding  series 
of  rotational  levels,  characterized  by  a rotational  quantum  number  J. 

In  Fig.  2.3-1  it  is  easily  seen  why  electronic  and  vibrational 
transitions  of  a molecule  cannot  be  observed  in  pure  form.  For  each 
electronic  transition  of  a molecule  the  vibrational  and  rotational 
energies  change  simultaneously.  The  entire  aggregate  of  electronic 
transitions  gives  rise  to  the  electronic-vibrational-rotational  spec- 
trum, whose  line  strengths  are  determined  By  the  energy-level  distri- 
butions of  the  molecules  and  by  the  proBabilities  of  the  corresponding 
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Figure  2.3-1. 


Diagram  of  Vibrational  and  Rotational  Levels  of  Two 
Electronic  States  of  a Molecule  [Adapted  from  [Orchin  ] 


13 


The  number  of  different  vibrational  transitions  accompanying  a 
change  in  electron  energy  and  the  number  of  different  rotational  tran- 


i 


i 


sitions  associated  with  a change  in  the  vibrational  energy  of  a 
molecule  may  be  very  large.  This  is  responsible  for  the  extremely 
complex  structure  of  molecular  spectra.  Additionally  vibrational 
or  rotational  levels  of  a molecule  may  overlap  with  each  other. 

Figure  2.3-2  shows  the  vibrational  modes  of  carbon  dioxide  and  water 
vapor,  two  of  the  atmospheric  consituents  that  have  absorption  lines 
in  the  Infrared  spectrum. 

Detailed  discussion  of  the  energy  level  selection  rules  may  be 
25  15 

found  in  Kroto  or  Goody  . The  interaction  of  electronic,  vibra- 
tional, and  rotational  energies  may,  as  a first  approximation,  be 
considered  as  additive.  Thus,  if  the  energies  were  strictly  inde- 
pendent, the  selection  rules  for  a combined  transition  could  be 
regarded  as  a combination  of  the  three  separate  selection  rules. 

Atoms  and  molecules  can  exhibit  electronic  line  spectra,  but  the 
overall  complexity,  overlapping,  and  simultaneous  changes  of  all 
three  forms  of  quantized  energy  of  molecules  requires  band  systems. 

Lastly,  we  consider  a fundamental  property  under  the  classical 
theory  of  a molecule.  That  is  the  electric  dipole  moment.  This 
moment  is  defined  as  the  product  of  the  magnitude  of  an  electric 
charge  and  the  distance  between  it  and  its  opposite  charge.  It  is 
this  property  by  which  absorption  occurs.  The  electromagnetic  radia- 
tion absorption  happens  when  an  electric  field  Interacts  with  the  motion 
of  the  molecule,  changing  the  dipole  moment  and  bringing  about  an 


14 


SYMMETRIC  STRETCHING 


NO  INFRARED  ACTIVITY  X>  2.738^in 


ANTI -SYMMETRIC  STRETCHING 


BENDING 


9 


9 


X s|4.986 


Xe  6.270 /xm 


18, 


Figure  2.3-2.  Vibrational  Modes  of  CO2  and  H2O  [Hudson  ] 
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energy  transition.  In  symmetric  molecules  some  vibrational  modes 
produce  no  change  in  dipole  moment.  As  a consequence,  no  Infrared 
activity  occurs. 


r 


2.4  Spectral  Line  Broadening 

Thus  far  in  the  discussion  of  spectral  lines  it  has  been  stated 
that  the  position  (frequency)  of  an  absorbing  line  due  to  an  allowed 
energy  transition  may  be  predicted  exactly  by  Eq.  2.3-2.  While  In 
theory  this  is  true,  in  practice  the  line  that  is  observed  is  not  a 
singled  frequency  slit  from  the  electromagnetic  spectrum.  The  line 
is  widened  by  a phenomenon  termed  line  broadening.  The  expansion  is 

4 

due  to  a variety  of  reasons  [Breene  ] , but  only  three  will  be  dis- 
cussed here.  Natural  line-broadening,  Doppler  line-broadening,  and 
pressure  or  Lorentz  line-broadening  are  the  main  mechanisms  considered 
for  atmospheric  calculations  and  are  therefore,  developed  below.  The 
line  breath  is  conventionally  taken  as  the  spectral  width  corresponding 
to  the  half  maximum  intensity  points  on  the  curve  describing  the  line. 
The  profiles  are  assumed  symmetric  hence,  half  the  line  breath  Is 
used  in  characterizing  the  line  shape,  and  Is  termed  the  half  width. 
Figure  2.4-1  depicts  the  Doppler  line  shape,  the  Lorentz  line  shape 
and  the  half  width  of  each. 

2.4.1  Natural  Broadening 

The  uncertainty  principle  states  that  the  product  of  the  uncer- 
tainty of  our  knowledge  of  the  energy  and  mean  life  of  an  atomic 
state  cannot  be  less  than  Planck's  constant  divided  by  4ir.  Since 
every  excited  state  of  a molecule  has  a finite  lifetime,  its  energy 
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f ( W - Wo) 
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cannot  be  measured  exactly.  Thus,  because  we  do  not  know  precisely 

a molecule's  energy  transition,  we  cannot  know  absolutely  the  f re- 

22 

quency  of  emmlssion  or  absorption  [Jamieson,  et  al  ] . This  uncer- 
tainty leads  to  a lower  limit  of  the  width  of  an  observed  spectral 
line,  termed  the  natural  width  (otjj)  • It  is  given  by 


X 

“n  ” 


(2.4-1) 


where  is  the  lifetime  of  an  energy  transition.  For  a typical 
vibrational  transition  is  of  the  order  of  0.1  seconds  [Goody^^]. 
Thus,  from  Eq.  2.4-1,  is  of  the  order  of  3*10  cm  in  wave 
numbers.  This  is  trivial  when  compared  to  other  line  widths  and  is 
neglected  for  atmospheric  calculations. 


2.4.2  Doppler  Broadening 

If  a monochromatic  radiation  source  is  moving  with  velocity 
whose  component  in  the  line  of  sight  to  the  observer  is  u,  the  fre- 
quency appears  shifted  by  an  amount 


Av  - - V 
c o 


(2.4-2) 


compared  with  the  frequency  of  the  source  when  it  is  at  rest (v^) . 
This  change  is  termed  the  Doppler  Effect. 

For  a stratified  atmosphere  it  may  be  assumed  that  the  transla- 
tional states  of  molecules  are  in  thermodynamic  equilibrium  and  their 
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velocity  u,  relative  to  an  observer.  Is  given  by  the  Maxwell  velocity 
26 

distribution  [Kuhn  ] 


p(u) 


{lirY  -p  (-  Ife  ) ■ 


(2.4-3) 


where 


m Is  the  mass  of  the  molecule, 
k Is  Boltzman's  constant, 
and  T Is  temperature  In  °K. 

When  u«c  this  probability  function  for  velocity  manifests 
itself  in  a broadening  of  the  spectral  line  observed. 

2 

Defining  the  doppler  half  width  at  half  maximum  [Anding  ] as 


D c 


^ ^2kTln2^^^^ 


,1/2 

= 3.58-10  1 V 


■’  (O' 


(2.4-4) 


where  M is  molecular  weight,  the  line  shape  is  predicted  by 

. . /ln2V''^  1 f . 


(2.4-5) 


15  —2  —1 
Goody  gives  the  two  extreme  values  for  as  3.3-10  cm  and 

-4  -1 

3.5-10  cm  . 
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2.4.3  Pressure  Broadening 


The  Interaction  of  a radiating  particle  with  another  particle 
has  an  effect  on  the  spectral  line  which  la  very  complex  In  nature. 

In  order  to  make  the  problem  manageable  it  is  necessary  to  Introduce 
simplifying  assumptions  whose  justifications  are  not  often  obvious 
and  whose  validity  is  limited.  The  theory  of  pressure  broadening 
has  developed  from  two  opposite  viewpoints.  Following  Lorentz,  it 
is  considered  that  a collision  breaks  the  otherwise  unperturbed  wave 
train  of  a radiating  particle.  The  line  width  Is  then  determined, 
not  by  the  natural  life  time  T,  but  by  the  shorter  time  T^  between 
collisions.  A fourier  analysis  of  the  wave  train  leads  to  a line 
width  proportional  to  the  effect  is  a kinetic  one,  depending 

on  the  relative  velocities  of  the  molecules. 

The  statistical  approach  was  first  applied  by  Holtsmark^^  to  the 
effect  of  ions  on  each  other.  It  was  later  applied  by  others  to  the 
mutual  effects  of  neutral  atoms.  It  neglects  the  motion  of  the  atoms 
and  regards  them  as  being  at  rest,  but  distributed  randomly  in  space. 
The  characteristic  frequency  of  the  radiating  or  absorbing  molecule 
is  changed  by  the  presence  of  another  molecule  according  to  some 
function  of  the  distance  separating  them.  Each  molecule  can  be 
regarded  as  emitting  one  frequency  at  any  given  time,  but  the  random 
effect  of  many  radiating  molecules  leads  to  a continuous  Intensity 
distribution  whose  calculation  is  purely  statistical  in  nature. 

The  Fourier  analysis  of  the  emitted  radiation  of  a particle  that 

has  been  interrupted  by  a collision,  leads  to  an  expression  for  the 

4 

intensity  of  the  disturbed  wave  function  [Breene  ] of 
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I(v) 


Cl  R J |”dt  exp 

r 


(2Tri(v  -v)t) 
o 


exp  2Tri  n^t) 


(2.4-6) 


where  I(v)  is  intensity  as  a function  of  v, 

R is  the  operation  of  taking  the  real  part  of  the 
expression, 

n(t)  is  the  phase  shift  due  to  interaction  as  a 
function  of  time,  and  the  bar  over  the  function 
indicates  a time  averaging. 


Now  if  there  is  negligible  phase  shift  before  time  t,  then  an 
encounter  occurs  and  the  integral  terminates,  we  have 


sin  2'rr(v  -v)t 
o 

2ir(v  -v) 
o 


(2.4-7) 


This  result  is  now  averaged  over  all  possible  values  of  t,  while 
neglecting  the  time  spent  in  collision.  The  interruption  theory  has 
now  become  Impact  theory. 


21 


According  to  the  kinetic  theory  of  gases,  the  distance  (1) 
traveled  between  collisions  by  a molecule  with  velocity  u follows 
the  distribution  function 


p(l)  dl  = — exp(-l/l  ) , 

1 “ 
u 


(2.4-8) 


where  1^  is  the  mean  free-path  for  molecules  of  velocity  u.  The 
distribution  function  for  the  time  between  collisions  follows  from 
the  substitution  udl  = dt,  and  is  given  by 


pCt)  dt 


(2.4-9) 


where  is  again  the  mean  time  between  collisions.  Now  the  product 
of  Eq.  2.4-9  and  Eq.  2.4-7,  integrated  over  all  time  yields  [Goody^^] 


I(v,u) 


which  if  normalized  so  that  its  area  for  all  v is  unity,  may  be 
written  as 


f(v,u)  « i 5 • (2.4-10) 

(v-v^)  + aj^(u) 


where  the  Lorentz  line  width  (a  ) is  defined  as 

L« 
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(2.4-11) 


2^T 


The  LorentE  line  width  is  very  strongly  dependent  upon  pressure 
54 

and  temperature  [Zuev  ] and  is  often  approximated  by 


(2.4-12) 


The  subscript  zero  denotes  some  reference  width  calculated  at 
the  reference  pressure  and  temperature  T^.  The  exponent  n is  an 
empirically  derived  number,  usually  set  to  0.5. 

Having  replaced  the  velocity  dependence  of  a with  a pressure, 

Li 

temperature  relation,  the  common  form  of  the  Lorentz  line  shape 
(Pressure  Broadening)  is  given  as 

^ 2 • 

(v-Vo)  + 

Typical  values  for  range  from  0.055  cm  ^ for  methane  to  0.11  cm  ^ 
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for  ozone  [McClatchey,  a^  ] . 

2.5  Summary 

A knowledge  of  the  composition  and  properties  of  the  atmosphere, 
together  with  an  understanding  of  molecular  energy  transitions,  and 
spectral  line  broadening  are  the  tools  to  be  used  for  atmospheric 
absorption  computations.  It  is  possible  to  identify  those  molecules 
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that  absorb  radiation  and  determine,  to  some  degree  of  accuracy,  the 

center  frequency  absorbed.  It  Is  also  possible  to  roughly  predict 

the  shape  of  the  spectral  line,  given  Information  on  the  composition 

of  the  atmosphere.  This  Information,  when  calculated  accurately,  may 

be  used  for  the  temperature  retrieval  scheme  alluded  to  in  the  Intro- 
33 

duct ion  [McClatchey  ] . 
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CHAPTER  III 

TRANSMITTANCE  AND  THE  VOIGT  FUNCTION 
IN  THE  FIFTEEN  MICROMETER  BAND 


3.1  Preface 

This  Chapter  addresses  the  subject  of  atmospheric  transmittance 
with  the  combined  effects  of  Doppler  and  Pressure  line  broadening. 

The  definitions  of  absorption  and  transmittance  are  first  presented 
then  a function  Including  both  types  of  spectral  line  spreading  la 
developed.  This  function  Is  then  Introduced  Into  the  transmittance 
differential  equation  and  the  final  form  of  the  transmittance  rela- 
tion is  given.  Next  the  frequency  band  used  here  is  briefly  discussed, 
together  with  the  validity  of  the  assumptions  made  that  limit  Inves- 
tigation to  this  Interval.  Finally  a summary  of  the  Important  aspects 
of  the  chapter  Is  given. 

3.2  Absorption  and  Transmittance 

The  interaction  of  radiation  with  matter  can  lead  to  three  types 
of  changes  In  the  energy.  (1)  Reflection,  in  which  a part  of  the 
Incident  energy  is  deflected  back.  The  reflected  radiation  as  a 
fraction  of  that  Incident  Is  the  power  of  reflection  of  the  substance. 
(2)  Scattering,  a part  of  the  radiation  energy  Incident  at  a definite 
direction  is  scattered  in  all  directions.  This  is  caused  by  an  In- 
teraction depending  on  the  relationship  of  the  radiation  wave  length 
to  the  scattering  particle  size  or  certain  optical  flaws,  such  as 
differences  in  refractive  Index  in  the  material.  (3)  Absorption,  by 
which  Is  designated  the  conversion  of  radiation  energy  Into  heat.  The 
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portion  of  the  energy  transformed,  expressed  as  a fraction  of  the 

incident  energy,  is  known  as  the  power  of  absorption.  Chapter  two 

developed  the  molecular  model  used  for  this  energy  absorption  and 

now  the  model  to  predict  absorption  not  on  a molecular  level,  but  a 

quantitative  level,  is  developed. 

These  three  effects  depend  in  a characteristic  manner  on  the 

wave  length  of  the  radiation  and  also  on  the  specific  properties  tf 

the  substance  to  be  investigated.  For  the  Infrared  region,  only  the 
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last  effect  is  important  [La  Rocca  and  Turner  ] . Therefore  it  will 
be  discussed  exclusively. 

The  radiant  energy  of  a definite  monochromatic  frequency  v 
passing  into  a body  of  length  Ax,  as  shown  in  Fig.  3.2-1,  is  denoted 
The  emanating  energy  is  reduced  by  some  amount,  proportional 
to  the  length  Ax,  and  is  denoted  I^(X2). 


For  atmospheric  purposes  it  is  recognized  that  the  reduction  in 
Intensity  of  the  radiation  is  due  to  the  absorption  of  a portion  of 
the  energy  by  the  molecules  present  In  the  path.  Figure  3.2-1  Is 
expressed  mathematically  as 


I^(X2)  = AX 


or  in  differential  form 


dl^(X)  = -I^(X)  adX. 


(3.2-1) 


The  constant  of  proportionality  is  termed  a and  is  also  called 
the  extinction  coefficient.  The  negative  sign  is  due  to  the  reduction 
in  Intensity.  Equation  3.2-1  is  integrated  to  obtain 


i^(*) 


I„(o) 


<1 


-adx. 


and  finally 


I (x) 

V 

I (o) 

V 


I^(x)  = I^(o)  exp (-ax) 


(3.2-2) 


Equation  3.2-2  is  known  as  the  absorption  law  and  is  usually 


named  after  Lambert-Bouguer . The  index  ax  is  called  the  optical 
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density  [Brugel  ] . The  optical  density  is  an  Index  of  a power  and 
must  be  dimensionless.  The  layer  thickness  x is  usually  measured 
in  centimeters,  thus  the  extinction  coefficient  is  expressed  in 


A graphical  representation  of  the  optical  density  (ax) , or  the 
radiation  energy  transmitted  (I^(x)),  both  for  a constant  layer 
thickness  (x) ; or  the  extinction  coefficient  (a)  as  a function  of 
frequency,  give  the  spectrum  of  the  substance  investigated.  However, 
other  quantities  are  usually  used  in  the  representation  of  spectra, 
namely  the  transmittance  (t)  or  the  absorption  (A) , both  given  as  a 
percentage  of  the  incident  radiation.  The  first  is  defined  as  the 
ratio  of  the  energy  passed  to  that  which  was  Incident.  So  from 


Eq.  3.2-2, 


I^(x) 

1^(0) 


= exp (-ax), 


(3.2-3) 


The  absorption  (A)  is  defined  as 


A = 1-T  = 1-exp (-ax). 


Accordingly,  the  extinction  (E)  is  defined 


(3.2-4) 


-Int  In 


1^(0) 

I^(x) 


(3.2-5) 


The  law  of  absorption  (Eq.  3.2-2)  gives  a first  quantitative  unit  for 
spectral  measurements.  Obviously,  the  absorption  of  a body  is  larger 
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the  greater  the  thickness  x,  or  the  larger  the  extinction  coefficient 
a,  for  the  wavelength  in  question.  In  the  last  resort,  however,  the 
number  of  molecules  absorbing  in  the  volume  is  what  is  responsible 
for  the  absorption.  Accordingly  the  extinction  coefficient  suffices 
completely  for  the  description  of  the  absorption  power  of  a sub- 
stance provided  the  number  of  molecules  in  the  volume  under  consid- 
eration remains  constant.  Unfortunately  this  does  not  apply  to  the 
atmosphere  where  the  number  of  molecules  per  unit  volume  is  dependent 
upon  the  pressure.  Therefore  a relation  must  be  derived  between  the 
concentration  of  a substance  and  the  absorption  it  causes.  Since  In 
the  atmosphere  the  number  of  molecules  responsible  for  the  absorp- 
tion depends  directly  on  the  pressure,  it  is  logical  to  write  the 
extinction  as  a linear  function  of  the  concentration  (or  density)  of 
the  absorbing  gas. 

E = K p(x)dx  (3.2-6) 

This  equation  is  referred  to  as  Beer’s  Law,  in  honor  of  its  discoverer, 
A.  Beer.  Now  from  Eqs.  3.2-6  and  3.2-5  we  obtain 

T = exp(-Kp (x)dx) . (3.2-7) 

This  combined  relationship  of  Eq.  3.2-6  and  Eq.  3.2-5  is  known  as  the 
Lambert-Beer  Law. 

For  mixtures  of  absorbing  materials  the  Incremental  absorption 
is  additive  in  the  absence  of  chemical  interaction  [Whif fen^^] . Thus 
Eqs.  3.2-2  and  3.2-7  may  be  rewritten  as 
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(3.2-8) 


1 (x)  •»  I (o)  exp(Z-a  x), 

j J 


and 


T = exp(-EK  p . (x)dx) , 

i ^ ^ 


(3.2-9) 


where  the  sumatlon  over  j accounts  for  the  different  gases  In  the 
atmosphere.  Equation  3.2-9  is  the  basis  of  the  whole  Infrared  spec- 
tral analysis.  It  is  true  in  most  cases  for  monochromatic  radiation. 
However  care  must  be  taken  when  it  is  applied  to  a band,  as  must  be 
done  with  any  measurement  instrument.  Brugel^  contains  a discussion 
of  the  limits  of  Eq.  3.2-9. 

3.2.1  The  Transmittance  Differential  Equation 

While  Eq.  3.2-9  is  correct,  it  is  not  in  a convenient  form  for 
calculation  of  atmospheric  transmittances  from  the  top  of  the  atmos- 
phere to  the  bottom  (ground  level).  Beginning  with  a form  of  Eq. 

3.2-1  and  utilizing  Beer's  Law  the  starting  relation  is  given  as 

dl  - -I  Kp(x)dx.  (3.2-10) 

The  expressed  dependence  of  intensity  (I)  and  the  absorption 
coefficient  (K)  on  wavelength,  temperature,  and  pressure  is  dropped, 
but  still  understood.  Rearranging  3.2-10  we  obtain 

dl  = -Kp (x)dx. 

I 

The  hydrostatic  equation  (2,2-1)  relates  density  (p)  to  pressure  (P) 
in  the  manner  desired.  It  is  restated  below  without  the  negative 
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sign.  For  atmospheric  considerations  we  begin  at  altitude  Z and 
pressure  P and  move  down  to  altitude  Z-dZ  and  pressure  P+AP.  Also 
the  mixing  ratio  M(P)  is  introduced  to  reflect  a relation  between 
total  pressure  and  the  pressure  of  only  one  gas  under  consideration. 
Chapter  five  contains  a detailed  explanation  of  Eq.  3.2--11. 
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g 


M(P)p(x)dx 


(3.2-11) 


This  yields 


da 
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m 

g 


dP, 


(3.2-12) 


where  the  dependence  of  the  mixing  ratio  on  pressure  is  understood. 

Instruments  measure  that  portion  of  the  intensity  of  radiation 
directed  downward  only  and  what  needs  to  be  considered  is  the  total 
intensity.  If  6 is  the  solar  zenith  angle,  from  Fig.  3.2-2 
Eq.  3.2-12  becomes 
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Figure  3.2-2.  Solar  Zenith  Angle  6 And 
Downward  Intensity 
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Integrating  this  equation  from  the  top  of  the  atmosphere  (P-0,  I-I  ) 

o 

to  a point  where  the  pressure  Is  P and  the  Intensity  I,  yields, 

I=I  P 

dl  sec  0 


g 


MKdP, 


I-I 


o 


sec  6 
g 


MKdP, 


or  from  the  definition  of  transmittance, 


T - exp  ^ 


sec  9 
g 


KMdP 


(3.2-13) 


Now  at  a reference  point  (T  , P ) calculate  a reference  transmittance 

o o 


T and  let 
o 


Thus, 


sec  9 
g 


- exp  < 


(3.2-14) 


-c 


KMdP  \ , 


and  at  any  other  position  given  by  P = P^+AP,  calculate 


32 


exp 


/'  -r 


AP 


N 

I 


Kt:Jv  ' 

J 


Dividing  this  expression  by  the  reference  transmittance  gives 


f 

I 

+AP  p n 

■N 

i 

0 

0 

) _c 

KMdp  - 

KMdP 

) 

1 

1 

, 

1 

0 0 

1 

V 

-j 

J 

P +AP 
o 


T / „ 

— = exp  \ -c 

o 


KMdP 


P 

V o 


Now  let  P “ P , a reference  pressure 
o r 


T ■ T , the  associated  reference  transmittance 
o r 


P 


T 


P +AP , and 
o 


P 

f 


T 

r 


exp  < -c 


KMdP  > . 


V 


(3.2-15) 


We  may  place  this  relation  In  (.  fferentlal  equation  form. 
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i 


Using  Leibnitz's  rule  for  differentiating  an  Integral 


1 ^ 
T dP 


-c 


p 

’ . . dP 


Since  K Is  a function  of  v,  the  last  two  terms  follow. 

dP 

nor  P are  functions  of  v,  so  — - ■ 0. 

r ’ dv  dv 


ISplegel^^], 


Neither  P 
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f 


p 


f • 


i ^ 

T dP 


-c 


5p  (KM) dp 


Using  the  fundamental  theorm  of  calculus  [Britten,  et  al^] , 


l.e.:  If  f (x)  = ^ g(x) 


then 


f (x)dx 


^ g(x)dx  - g(x). 


subject  to  boundary  conditions; 


- - -cKM. 

T dp 


Writing  this  expression  to  explicitly  show  the  functional  dependence 
upon  frequency,  temperature,  and  pressure  and  utilizing  the  knowledge 
that  the  extinction  coefficient  (K)  of  each  gas  may  be  superimposed 
(as  In  Eq.  3.2-9),  the  transmittance  differential  equation  Is  obtained. 


dT(v.T.P) 

dP 


sec  6 
g 


^Kj(v,T,P)Mj(P)^  t(v,T,P).  (3.2-16) 


The  Iterative  solution  of  this  expression  yields  the  desired 
transmittance  values  from  the  top  of  the  atmosphere  to  ground  level. 
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w 


3.2.2  The  Absorption  Coefficient  K(v.T.P) 

For  a given  gas  there  exists  a variety  of  forms  for  the 

30 

absorption  coefficient  K(v,T,P)  [L'Vov  ].  In  general  the 
coefficient  can  be  written 


Cj(v,T,P) 


S f(v-v  ). 


(3.2-17) 


S.  represents  the  strength  of  the  1th  spectral  line,  v Is  the 
1 o^ 

center  frequency  of  the  line  and  the  function  f(v-v^  ) defines  the 
line  shape  of  the  1th  line.  In  theory  the  contribution  of  each  spec- 
tral line  Is  added  together  (the  summation  over  1)  to  obtain  an 
absorption  coefficient  for  a molecular  species  (gas  type).  Several 
authors  give  expressions  for  the  line  Intensity  (S^)  [Goody^^, 
Hertzberg^^] . Miller,  e£  al.^^  gives 


S Q (T) 
IJ  J 


T 


(3.2-18) 


where  S Is  a reference  strength  calculated  at  (T  ,P  ) , 

o o o 

E Is  the  energy  of  the  lower  state  Involved  In  the 
molecular  transition, 
k Is  the  Boltzman's  constant, 

T Is  temperature  in  degrees  Kelvin,  and 
Q Is  the  vibrational  partition  function. 

This  strength  function  Is  one  that  currently  being  used  for  atmospheric 


transmittance  calculations  and  Is  adopted  here. 


36 


What  remains  of  Eq.  3.2-17  to  be  specified  is  the  line  shape 

function.  Chapter  two  enumerated  the  three  most  important  line 

broadening  mechanisms.  In  practice  the  investigator  must  decide  if 

the  problem  at  hand  can  be  resolved  by  a single  broadening  function 

or  if  a more  complicated  one  is  needed  to  accurately  describe  the 

absorption  coefficient.  There  are  several  approaches  to  the  problem 

35 

and  the  interested  individual  is  referred  to  Miller,  et  , 
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Andlng  , or  La  Rocca  and  Turner  . For  rapid  line-by-line  computation 
of  transmittance  at  the  sacrifice  of  some  accuracy  the  Lorentz  line 


function  employed  in  this  work  is  developed  in  the  next  section. 

3.3  The  Voigt  Function 

The  first  twenty  years  of  the  twentieth  century  were  marked  by 

considerable  achievements  in  the  field  of  atomic  absorption  theory. 

During  this  period  the  fundamental  relationships  between  absorption 

and  the  atomic  constants  were  established  and,  as  was  discussed  In 

Chapter  two,  the  theory  for  pressure  broadening  of  lines  was  formu- 
29 

lated  [Lorentz  ] . The  theory  of  Natural  and  Doppler  broadening  had 

been  already  developed.  While  several  individuals  realized  that  line 

broadening  phenomena  could  not  be  completely  predicted  unless  all 

48 

causes  were  accounted  for,  it  was  W.  Voigt  who  first  proposed  that 
the  effects  of  Doppler  and  Collision  broadening  could  be  dealt  with 
in  a single  line  shape  function.  Today  this  relation  bears  his  name. 

The  relation  may  be  derived  by  assuming  that  the  line  shape 
function  follows  a Doppler  profile  near  the  center  and  Lorentz  in 
the  wings. 
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The  Lorentz  line  shape  was  developed  in  Chapter  two  and  is 
restated  here. 


1 L 

f(v-v  ) = i ^ 

(v-v  ) +a 

O L 


(3.3-1) 


The  Doppler  effect  that  causes  a frequency  shift  of  the  monochromatic 
wavelength  » emitted  by  a molecule  at  rest,  is  written 


V'  = v_^(l  - P 

O O c 


(3.3-2) 


Again  u is  the  line-of-slght  velocity  component  of  the  molecule. 
From  Eq,  3.3-2  two  relations  to  be  used  later  are  obtained. 


C^(v'  - V )^ 
2 o o 


(3.3-3) 


du  = — dv 
V o 
o 


(3.3-4) 


Now  if  the  translational  states  of  the  molecules  under  consideration 
are  in  thermodyriamld  equilibrium, the  fractional  number  of  molecules 
with  velocities  between  u and  u + du  is  given  by  the  Maxwell  velocity 
distribution  p(u). 


(3.3-5) 


Substituting  Eq.  3.3-3  and  3.3-4  into  3.3-5  yields 


p(u)du 


2/  - •>2 
me  (v^  - v) 
o o 


C , , 
— dv  . 
V o 
o 


(3.3-6) 
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The  contribution  to  the  shaping  function  from  all  Doppler  shifted 
components  will  be  the  probability  function  p(u)du  times  the  Lorentz 
line  shape.  Integrated  over  all  frequencies.  It  may  be  expressed  as 


f(v  - v^) 


p(u)  f'(v  - V )du. 

o 


(3.3-7) 


The  function  f'(v  - v^)  Is  given  by  Eq.  3.3-1  with  the  center 
frequency  of  absorption  (v^)  replaced  by  the  Doppler  shifted  frequency 


(vp.  Thus, 


and 


f'Cv  - V ) » — 
o 7r 


2 2 * 
(V  - + of 

o L 


(3.3-8) 


V-V  =(v-v)  - (\)'  - V ). 
o o o o' 


(3.3-9) 


Now  utilizing  Eqs.  3.3-9,  3.3-8,  3.3-6  and  3.3-7  the  combined 
Doppler-Lorentz  function  Is  obtained  [Volgt^^]. 


f(v-v  ) - 
o 


1 

IT 


“ — 

o^  exp  - 

2 - 2 
me  (v  -V  ) 

0 0 

2kTv  ^ 

L o J 

[(V-V^) 


iTTkTj 


c_,^. 


Defining  the  variable 


'’o  « 


(3.3-10) 


2^ 

1/2 

fv'-v  \ 

me 

o o 1 

2kT 

"o 

i 

1 ° / 

(3.3-11) 
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we  have  immediately , 


and 


(3.3-12) 


(3.3-13) 


So  now  Eq.  3.3-10  becomes 


f(v-v^) 


(3.3-14) 


exp (-t  ) dt 


X f. 


.4  1/2 


j ) / r (v-v  ) - ( V "-v  ) ] ' 

i 2kT  I S ^ ^ o'^  ^ o 


(3.3-15) 


Now  using  Eq . 3.3—11  again  and  defining 


(3.3-16) 


(3.3-17) 

(3.3-18) 


where  cx^  is 


the  Doppler  half  width,  the  relation  becomes 


f(v-v^) 


2 

y exp(-t  )dt 

2 2 * 
[x-tr  + y 


(3.3-19) 


The  absorption  coefficient  for  the  mixed  Doppler-Lorentz  (Voigt) 
function  is  then  given  by  Eq.  3.2-17  and  the  following  new  defini- 
tion of  strength 


(3.3-20) 


where  is  defined  in  Eq. 


3.2-18. 


Thus 


Kj(v,T.P) 


2 

y exp(-t  )dt 

2 2 
y + (x-t)^ 


(3.3-21) 
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Other  derivations  of  this  Voigt  profile  are  given  in  Penner  , Goody 
or  more  recently  Young^^. 


Combining  Eq.  3.3-21  and  the  transmittance  differential 
equation  (3.2-16)  one  obtains  the  relation  for  atmospheric  transmit- 
tance considering  the  Voigt  profile.  Several  authors  assert  that 

this  function  should  be  used  if  completely  accurate  results  are  to  be 
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obtained  {LaKocca  and  Turner  , Andlng  , and  Rodgers  ] . The  sum- 
mary of  this  chapter  collects  all  the  equations  thus  far  derived 
that  are  needed.  The  evaluation  of  the  Voigt  profile  Is  discussed 
In  the  next  chapter. 

Figure  3.3-1  shows  the  relation  between  the  Doppler  line  shape, 
Lorentz  line  shape  and  the  Voigt  profile. 

3.4  The  Fifteen  Micron  Band 

Figure  3.4-1  shows  the  electromagnetic  spectrum  and  the 
nomenclature  used  for  different  regions.  This  work  Is  limited  to 
the  so  called  fifteen  micron  band.  The  limits  of  the  band  are  not 
well  defined  and  in  fact  depend  much  on  the  instrument  used  to 
Investigate  the  region.  As  the  name  implies  however  It  Is  generally 
around  the  ISpm  portion  of  the  spectrum. 

The  NIMBUS  6 High  Resolution  Infrared  Radiation  Sounder  (HIRS) 

Is  a third  generation  Infrared  radiation  sounding  satellite  possessing 
many  new  features  for  greatly  Improved  capability  of  sounding  the 
earth's  atmosphere.  The  Investigations  of  this  work  used  the  speci- 
fications of  the  NIMBUS  6 to  calculate  the  transmit tances  needed  to 
provide  correct  temperature  soundings.  NIMBUS  6 uses  seven  channels 
in  the  fifteen  micron  band.  The  channels  are  monitored  by  a multi- 
channel filter  radiometer.  The  band  provides  better  sensitivity  to 
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STRENGTH 


Figure  3.3-1.  Doppler,  Lorentz  and  Voigt  line  shapes. 
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I 

i 


the  teaperature  of  the  relatively  cold  regions  of  the  ataosphere 
43 

[Slssala  ] than  can  be  achieved  with  the  4.3ya  band  channels. 

These  cooler  temperatures  correspond  to  higher  altitudes  and  lover 
pressures.  At  these  levels  the  Lorentz  line  width  Is  reduced  and 
the  Doppler  line  width  Is  broadened.  Therefore  for  most  accurate 
transmittance  calculations  both  line  broadening  schemes  must  be 
represented.  For  this  reason  the  Voigt  profile  Is  used  In  calcula- 
tions . 

The  Important  engineering  aspects  of  NIMBUS  6 design  are 
depicted  In  Table  3.4-1  and  Fig.  3.4-2. 


Central 

Interval  Between 

Level  of 

Channel 

Wavenumber 

50%  Response  Points 

Primary 

Peak  Energy 

lhad>er 



^SS..  ? 

Absorber 

Contribution 

1 

668 

2.8 

CO2 

30mb 

2 

679 

13.7 

CO2 

60mb 

3 

690 

12.6 

CO2 

lOOab 

4 

702 

15.9 

CO2 

250nd> 

5 

716 

17.5 

n 

0 

10 

500nd> 

6 

733 

17.6 

CO2/H2O 

750Bd) 

7 

749 

18.4 

CO2/H2O 

900mb 

Table  3.4-1. 

NIMBUS  6 Characteristics, 

43 

(Adapted  from  Slssala  ) 

Figure  3.4-2  depicts  Che  filter  responses  for  the  seven  channels 
of  interest.  For  atmospheric  calculations  It  Is  necessary  to  consider 
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Flaure  3.4-2.  HIRS  filter  responses  [Sissala  ]. 


all  lines  whose  center  frequency  (v^)  is  within  a channel,  for  the 
transmittance  of  that  channel.  Thus,  for  example,  the  summation 
over  1 In  Eq.  3.3-21  runs  over  those  lines  whose  center  frequency 
Is  between  665.2  cm~^  and  670.8  cm  ^ for  channel  one.  In  addition 
It  may  be  necessary  to  consider  the  wings  of  lines  whose  center 
frequency  lies  outside  the  channel  Interval.  Figure  3.4-3  depicts 
such  a line  (wavenumber  ) . Although  the  center  frequency  Is  out- 
side the  channel  limits,  the  right  wing  of  the  line  contributes  to 
absorption  within  the  channel. 


3.5  Summary 

The  transmittance  differential  equation  to  be  used  In  this 
Investigation  Is  given  by 


d t(v.T.P) 
dP 


- 


sec  9 

g 


^Kj(v,T,P)  Mj(P) 

J 


T(v,T,P), 


(3.2-16) 


Kj  Is  called  the  absorption  coefficient  and  has  different 
representations,  depending  upon  the  nature  of  the  Investigation. 

For  Instances  where  the  combined  effects  of  Doppler-Lorentz  (Voigt) 
broadening  are  to  be  considered,  the  absorption  coefficient  Is 
expressed  by  the  Voigt  profile. 


2 

y exp(-t  )dt 
y^  + (x-t)^ 


(3.3-21) 


where 


( 

' ) 
1 

I < 
1 ; 
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la 

“d 


Is  the  line  strength. 


(3.3-20) 


(3.3-17) 

(3.3-18) 

(2.4-12) 

(2.4-4) 


NOTE; 


1/2 

The  doppler  width  o^dnl)  la  used  Instead  of  In  all 

2 

equations  lAndlng  ] . 


M Is  the  molecular  weight,  the  Lorentz  half  width,  the 

Doppler  half  width,  Q the  vibrational  partition  function  (see  Chapter 

six),  E the  energy  of  the  lower  state  of  a transition,  and  k Boltzaan's 

constant.  The  subscript  zero  Indicates  a reference  value  calculated 

at  (T  , P ),  except  when  lised  with  v.  Then  It  Indicates  the  center 
o o 

frequency  of  absorption.  The  subscript  j la  for  gas  type  and  the 
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subscript  1 is  for  the  absorbing  lines.  J runs  for  all  molecular 
species  of  concern,  while  the  1 counter  runs  for  all  the  lines  that 
fall  within  a channel  of  the  Instrument  used  for  measurements.  The 


T 


limits  of  the  channel  may  have  to  be  extended  by  some  amount  to 
account  for  the  wings  of  lines  whose  center  frequency  Is  outside 
the  channel.  The  channels  used  here  and  their  filter  responses  are 
depicted  In  Table  3.4-1  and  Fig.  3.4-2  respectively. 
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CHAFTER  ly 


CQMFUTER  ALGORITHMS  FOR  THE  VOIGT  PROFILE 

4.1  Preface 

The  Voigt  profile  developed  in  chapter  three  is  widely  accepted 
hut a until  recently,  not  widely  used  because  It  is  difficult  to 
qulcRly  evaluate.  This  chapter  presents  three  methods  of  evaluation 
that  are  easily  computerized.  The  first  was  developed  by  John  Kellkopf, 
the  second  by  S.  R.  Or ay son  and  the  last  is  a new  method  developed 
here.  In  searching  for  methods  to  compute  the  profile,  speed  and 
accuracy  were  the  primary  concerns.  In  a line-by-line  transmittance 
model  as  much  as  75%  of  the  computlonal  time  may  be  used  to  evaluate 
the  Voigt  function.  Any  reduction  in  calculation  time  for  the  function 
would  then  obviously  have  a dramatic  effect  on  the  transmittance  pro- 
gram execution  time. 

4 

The  accuracy  required  of  any  algorithm  Is  about  one  part  in  10  . 

The  uncertainty  of  absorption  line  parameters  limit  transmittance 
accuracy  so  that  this  degree  of  exactness  is  all  that  is  needed. 
Function  subprograms  for  each  method  are  listed  in  appendices. 

The  results  of  tests  for  speed  and  accuracy  of  each  method  are 
presented  last.  The  design  of  the  examination  is  specified  and  the 
outcome  is  tabulated.  The  method  selected  for  use  in  these  trans- 
mittance calculations  is  specified  and  some  comments  on  its  applica- 
bility are  given. 

4.2  J.  F.  Kielkopf's  Voigt  Algorithm 

The  Voigt  profile  derived  in  chapter  three  is  written  in  slightly 
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different  form  by  Klelkopf  yollowtng  the  development  of  hla 
algorithm  the  tvo  forms  vlll  be  related.  Starting  as  Klellcopf,  the 
Voigt  profile  ffla7  be  written 


expl-CvVB^)  ]dv" 

IsJ  + (v-vO^J 


(4.2-1) 


A standardised  form  of  the  function  that  does  not  change  Its 
width  or  height  as  a function  of  parameters  la  given  by 


V(B,,  eX) 

■ 1(6..  8.)  > 


(4.2-2) 


g 


where  v ■ gx  and  the  peak  value  of  V Is  1.  For  large  values  of  x 
the  function  may  be  approximated  by 


1+g 


(4.2-3) 


where 


1 ■ 3^/e  , 


g - gg/e  . 


(4.2-4) 

(4.2-5) 


The  letters  g and  I represent  the  gatisslan  and  lorentslan  parts 
of  the  line  shape  respectively.  Near  the  line  core,  gausslan  proper- 
ties appear  In  U and  an  approximation  to  the  standardised  function  Is 


nhAxe 


GCx)  " expl-(,iln2)x^]  , C4.2-7) 

and 

L(x)  - (4.2-8) 

1-ht^ 

A function  of  the  form  of  Eq.  4.2-6  has  substantial  errors  In  the 
Halt  as  X grovs  large,  however.  Therefore  corrections  and  adjust- 
ments are  made  to  the  relations  used  for  1,  g,  and  U. 

The  graphs  given  by  Van  de  Hulst  and  Reesnlck^^  show  that  the 
gausslan  and  lorentzlan  portions  of  a line  are  related  approximately 
by 

g^  - (l/ln2)(l-l).  (4.2-9) 

A small  quadratic  correction  results  in  a highly  accurate  relation. 

g^  - (l/ln2)[l-(l+eln2))l+(eln2))l^]  (4.2-10) 

e - 0.0990 

The  lorentzlan  fraction  Is  similarly 


i 


2 

Il+eln2+I  (l-cln2)^+(41n2/a^) 


Hera  a Is  called  the  Voigt  parameter  and  Is  defined  as 


a 


e 

g 


(4.2-11) 


(4.2-12) 


53 


Frcao  Eq,  4.2-4  we  have 


6 " 


(4.2-13) 


Next  what  la  needed  is  the  peak  Intensity  I.  It  Is  given  In  terns  of 
the  complementary  error  function  as 


1(0 


Z* 


exp (a  ) erfc(a). 


(4.2-14) 


There  are  a large  number  of  methods  to  compute  the  coiiq>lementary  error 
function.  The  one  chosen  here  Is 


erfc(a)  " f(a)  exp  (-a^) 


(4.2-15) 


where 


f (a) 


a>lt  + bjt^  + b3t^) 


l+bpa  • 

- 0.61658 
63  - 1.32554 


bp  - 0.47047 
bj  - -0.16994 


(4.2-16) 


(4.2-17) 


Confining  Eqs.  4.2-14,  15  and  16, 


U6,,  Bg) 


bit  + b2t 


(4.2-18) 
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Bilgismtag  vlth.  Eq.  4.2-6  and  adding  a correction  term,  the  final  form 

24 

of  UCn.:n)  wae  determined  bjr  Klelkopf  as 

uCnsx)  ■ (l-ri)GGt)+nLOc)+n(l-n)E(x)lG(x)-LO»)].  (4.2-19) 

The  function  E(x)  was  determined  empirically  as 


E(x) 


0.8029  - 0.4207x^ 

1 + 0.2030x^  + 0.07335x^ 


(4.2-20) 


Now  the  relation  between  the  variables  of  Eq.  4.2-1  and  the  standard 
form  of  the  Voigt  profile  Is  shown. 


f(v-v^) 


2 

y 


j 

+ (x-t) 


2 


dt. 


(4.2-21) 


where 


v-v 


“d 


° (ln2)^/^. 


(4.2-22) 


(4.2-23) 


Substituting  these  definitions  In  Eq.  4.2-21  yields 
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V 


Nov  let 


«.oa 


2 


+ I Cv-V  ) - 

o 


t - v76g  . 

(4.2-24) 

V-V  = V , 

0 ’ 

(4.2-25) 

(1.2)^'^  ■ 

(4.2-26) 

“l  “ * 

(4.2-27) 

so  that 

f(v-v  ) " — 

O IT 


f(v-V^) 


exp 

[■sr] 

dv' 

4 

+ (v-v')^ 

(4.2-28) 
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This  Is  exactly  the  form  given  b/  Klelkopf  , so  the  relations  of  Eqs. 


4.2-24,  25,  26,  and  27  are  those  desired. 

Nov  in  terms  of  the  line  parameters  prevlousl7  defined,  the  pro- 


1 


C4.2-29) 

(4.2-30) 

(4.2-31) 

(4.2-32) 

(4.2-33) 

(4.2-34) 


(4.2-35) 


Ci)  LOc) 


1 


(A. 2-36) 


CIO)  ECx) 


0.8029  - 0.A207X* 

1 + 0.203x^  + 0.07335x^ 


(4.2-37) 


(11)  u(n:x)  - (l-Ti)G(x)  + nL(x)  + tiC1-ti)E(x)Ig(x)  - LCK)1  (4,2-38) 

(12)  V(8^  e^.v)  - 0(n:x)  I(8ji,3g)  (4.2-39) 


Th«s«  tv*lv«  Btapa  yl«ld  the  d«alr«d  valua  of  tha  Voigt  proflla.  Tha 
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accuracy  of  thla  method  la  given  [Klelkopf  ] aa  10  tinea  the  peak 
intenaity.  The  TORTltAN  IV  function  aubprogram  implementing  thla 
method  la  Hated  in  Appendix  A. 

4.3  S.  R.  Drayaon’a  Voigt  Algorithm 

In  developing  a computer  algorithm  to  evaluate  the  Voigt  function, 
Drayaon  required  accuracy  of  about  one  part  in  10^.  Beyond  thla  the 
chief  aim  vaa  rapid  computation  of  tha  function,  ehich  ia  repeated 
here  for  clear ity. 


f(v-v^) 


r + (*-t)' 


(4.3-1) 


It  ta  •aatly'  acBn  from  the  definition  of  x that  It  can  aaeuma  values 
that  are  both,  positive  and  negative.  However  since  It  Is  assumed 
that  fCvov^}  Is  an  even  function  In  x,  l.e.  fCs)  " only  posi- 

tive values  of  x are  needed  for  line  broadening  problems.  Similarly 
It  Is  seen  that  y will  assume  only  positive  values  as  long  as  the 

Lorants  half-width,  (o^)  Is  Included.  Vor  other  applications  nega- 

13 

tlve  values  of  y occur  iFrled  and  Conte  ] . Limiting  x and  y so  that 
they  are  greater  than  or  equal  to  sero  defines  the  first  quadrant  of 
the  xy  plane. 

Hov  consider  the  complex  probability  function  defined  as 


,») . i 


(4.3-2) 


where  Z ■ x -f  iy,  x >.  o,  y > o.  Separating  this  function  into  real 
and  Imaginary  parts  yields 


w(Z)  - ± 


le  ^ dt 
(xfly)  - t 


1 

IT 




(x-t)  + ly  (x-t)  - ly 


1 

IT 


2 2 
V e~^  + l(x-t)e~^ 

(x-t)^  + y^ 


dt 
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-t 


Cx^tr  + y" 


dt 


+ - 
II 


M 


-t 


Ot-t)  + y 


dt 


C4.3-3) 


It  is  thus  shovn  that  the  Voigt  profile  Is  the  real  part  of  the  cob- 

plex  probability  function  v(Z}.  Several  authors  have  tabulated  values 

of  this  function  for  both  the  real  and  Imaginary  parts,  but  the  aost 

cooplete  account  of  Its  mathematical  properties  appears  to  be  that  of 

12 

Feddeyeva  and  Terantdv  or  more  recently,  for  atmospheric  applications, 

3 

Armstrong  . 

In  order  to  evaluate  the  real  part  of  w(Z)  as  rapidly  as  possible, 
the  xy  plane  is  partitioned  as  shown  In  Figure  4.3-1  [Drayson^^]  and 
different  approximating  functions  are  used  In  each  region. 

In  Region  I w(Z)  may  be  rewritten  as 


w(Z) 


F(Z) 


F(Z) 


(4.3-4) 


For  real  Z,  the  function  F(Z)  Is  known  as  Dawson's  function.  The 
coefficients  (d^)  of  a Taylor  series  expansion  of  F(Z)  about  any  real 
point  (x,o) , have  the  relationships 

‘‘o  " **1  " ^"^*‘*o*  Vl  • ■ ^ ‘‘n-l^ 

for  n ■ 1,  2,  3,  ...  iDawson^^). 
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r - ^ - ■ 

Thu«,  for  Region  I,  Davson'a  function  is  first  evaluated  on  the 

e-axls  and  a Tailor  series  is  used  to  find  the  litaglnary'  part  of  ?(Z) 

' and  thus  the  real  part  of  vCZ) . 

I Davson'a  function  may'  he  evaluated  by  series  expansion,  asyiqttotlc 

I 19 

\ expansion,  or  a Chehyshev  expansion  IHixmmer  ].  In  the  method  deve- 

loped for  his  program  Drayson  evaluates  Dawson's  function  Initially 
by  Chehyshev  expansion  on  the  x-axis  for  x ■ 0,1,  0.3,  0.5,  ...,  4.9. 

The  coefficients  d^,  62,  d^,  are  computed  at  each  point  using  Eq. 
4.3-5  and  stored  for  later  use.  This  process  is  required  only  once. 
Dawson's  function  is  then  evaluated  for  arbitrary  x by  a Taylor  series 
expansion  about  the  nearest  tabulated  point  (x  0.1,  0.3,  ...,  4.9). 

The  number  of  terms  needed  in  the  series  expansion  of  F(Z)  is  s function 
of  X and  y.  To  avoid  a test  for  convergence  of  the  series,  an  enpiricel 

function  was  developed  to  determine  the  number  of  terms  required. 

9 

For  Region  II  first  note  that  it  has  been  shown  by  Chehyshev  and 
39 

later  K.  A.  Posse  that  an  Integral  of  the  type 
fb 

dt 

Z-t 

' a 

can  be  represented  by  the  continued  fraction 


62 


f 


vtlftte  % and  o(  ai'e  ^ a recurrence  fcrmula.  Fried  and 

n n 

13 

Conte  applied  this  technique  to  the  function  v(Z)  to  obtain  Its 
fractional  expansion. 

1 1 

wCZ)  « 2 + 1/2 

Z + 1 

Z 3/2 

Z+... 


w(Z) 


1 . 11/2  1 3/2  ...  n/2  ... 

“T72  Z+Z  + Z+  Z Z+ 


(4.3>6) 


Along  the  line  y 1,  which  separates  Region  I and  II,  the  Taylor 
series  and  the  continued  fraction  require  about  the  same  computation- 
al time.  The  number  of  terms  needed  in  the  fractional  expansion  varies 
from  n <■  4 to  n « 19  [Drayson^^].  Again  an  empirical  function  was 
developed  to  determine  the  number  of  terms  required  and  thus  avoid  s 
check  for  convergence. 

For  Region  HI  the  complex  argument  (Z)  has  a sufficiently  large 
2 2 

modulus,  i.e.  |z|  - (x  -fy  ) , so  that  Gauss-Hermlte  quadrature  Is 

12 

effective  [Faddeyeva  and  Terentdv  ] . The  Voigt  profile  Is  approxi- 
mated by  a four-point  quadrature  In  Region  Ilia  and  a two— point  quadra- 
ture In  Region  Illb.  The  expressions  used  are 


Z 

IT 


2 

y 


e 

+ (s-t) 


2 


dt 


W 


^ y + (x-t^)" 


(4.3-7) 


In  Region  Ilia,  and 


1 

H 


1 


A 
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1 


(4.3-8) 


Z 

ir 


( 


2 

y 


dt  * 


In  Region  Illb.  Here  are  the  weight  factors  and  the  t^'a  are  the 

g 

roots  of  the  Hermlte  polynomial  of  degree  nfl.  Carnahan,  et  al. 
has  an  explanation  of  Gaussian  Quadrature  methods  as  well  as  tables 
of  and  for  the  Gauss-Hermlte  Quadrature  method. 

Utilising  the  four  regions  of  Fig.  4.3-1  and  the  three  methods 
described  to  approximate  the  real  part  of  Eq.  4.3-2,  Drayson  has 
written  a FORTRAN  IV  function  subprogram  that  he  claims  Is  accurate 

4 

to  about  one  part  In  10  and  much  more  rapid  than  any  other  currently 
available  to  evaluate  the  Voigt  profile.  A listing  of  the  progrsn  Is 
contained  In  Appendix  B. 


4.4  Proposed  Voigt  Algorithm 

This  method  begins,  as  does  Drayson* s with  the  coo^lex  probability 
function 


w(Z) 


(4.4-1) 


Again,  Z ■ X -f  ly  and  x,  y ^ o.  It  was  shown  In  section  4.3  that  the 
real  part  of  this  function  Is  the  Voigt  profile.  From  the  definition 
of  the  complementary  error  function  of  complex  argument. 


(4.4-2) 
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Thftrsfore 


vCZ)  ■ erfcC-iZ), 


or  In  terms  of  the  error  function  of  complex  argument 


C4.4-3) 


v(Z)  ■ e*^  |l  - erfC-lZ)  ) . 


(4.^4) 


It  is  possible  to  approximate  w(Z)  by  utilizing  Euler's  formula 
for  e and  a series  expansion  of  the  error  function  for  co^lex  ergu- 
sMnt.  Euler's  formula  applied  here  is 
2 

e ^ - exp(-x^  + y^)|co8(2xy)  - 1 sln(2xy)|  . (4.4-5) 
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The  series  expansion  for  the  error  function  is  [Salser  ] 


TT 

n-o 

Condtlning  Eq.  4.4-5  and  4.4-6  into  Eq.  4.4-4  yields  the  algorithm 
desired. 

It  was  noted,  however,  that  for  large  modulus  the  series  expan- 
sion of  the  error  function  required  upwards  of  50  terms  for  convergence. 
Therefore,  the  xy  plane  was  divided  in  a manner  similar  to  that  of 
Dreyson.  The  difference  is  that  there  are  only  three  regions  and  two 
approximating  functions.  Figure  4.4-1  depicts  the  three  regions. 

The  series  is  used  in  Region  I.  In  order  to  avoid  testing  for 
convergence  a computer  routine  was  developed  to  determine  the  number 
of  terms  of  the  series  required  for  various  values  of  x and  y in  the 


rtglon.  The  results  of  that  routine  are  shown  In  Fig.  4.4-2.  It 
vaa  noticed  that  if  the  points  of  x > 0.1  and  x - 2.0  are  used  (a 
rough  linear  fit),  the  lines  for  different  values  of  7 have  the  same 
slope.  The  line  Just  moves  up  for  Increasing  values  of  y.  Therefore 
tha  relation  used  to  fix  the  number  of  terms  was  chosen  to  be 

N - 6.842  X + 8.  (4.4-7) 

The  position  of  this  line  with  respect  to  the  test  results  Is  also 
shown  In  Fig.  4.2-2.  In  Implementation  the  value  of  N Is  truncated 
to  yield  an  Integer.  The  minimum  for  N was  set  to  8 to  assure  con- 
vergence. The  single  exception  la  for  x ■ o.  Then  the  number  of 
terms  Is  set  to  15. 

The  limits  of  the  regions  depicted  In  Fig.  4.4-1  were  chosen 
after  analysing  different  approximations  for  w(Z) . It  was  decided  that 
a quadrature  of  the  Integral  would  be  the  simplest  and  most  accurate 
for  regions  two  and  three.  Using  a modified  Gauss-Laguerre  quadra- 
ture, Abramowltz^  obtains  two  approximations  that  are  used  for  w(Z) . 

The  limits  expressed  by  Abramowltz  were  modified  since  a reduced  accu- 
racy from  that  which  he  specified  was  acceptable.  For  Region  II  a 
three  point  formula  yields 

w(Z)  - IZ  (-/-  ■ + + ■/—  ]»  (4.4-8) 

\Z-A2  Z-Afi/ 

- 0.4613135  A^  - 0.1901635 

Aj  - 0.9999216  A^  - 1.7844927 
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NUMBER 


- 0.002883894 


Ag  - 5.5253437 


Only  the  real  part  of  equation  4.4-8  la  required,  so  the  relation  la 
separated  Into  real  and  Imaginary  parts. 


wGS)  " iz 


+ 12xy  Ot^-y^-A^)  + 12xy  U^-y^-Ag)  + 12jty 


Aj^(x^-y^-A2)  - i2Aj^xy  A2Cx^-y^-A^)-12A2xy 
" ^ ^ ix^-Y^-k^)^  + (2xy)^  ^ Cx^-y2_A^)2  + (2xy)2 


Aj(x^-y^-Ag)  12AgXy 

(x^-y^-Ag)^  + (2xy)^ 


I j 


- IZ 


Aj^Cx^-y^-A^) 


A3(x^-y^-A^) 


AgCx^-y^-Ag) 


(x^-y^-A2)^+(2xy)^  (x^-y^~A^) V(2xy)^  (x^-y^-Ag)^+<2xy)^ 


R 


-1 


1Z(R-1I) 
(-yflx) (R-11) 
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Th»  real  part  of  this  relation  Is 


[ 

r 

i 

I 


r(w(Z)3  i -yR  + xI. 


In  terms  of  the  constants  and  x and  y 


R^v(Z)J  - 


o 2 ,2  2 . . 

2x  y-(x  -y  -A2)y 
(x^-y^-A2)^  + (2xy)^ 


+ A 


2x^y- (x^-y^-A^)y 
^ (x^-y^-A.)^  + (2xy)^ 


« 2 ,2  2 ^ . 

2x  y-(x  -y  -Ag)y 

(x^-y^-Ag)^  + (2xy)^  * 


(4.A-9) 


Equation  4.4-9  Is  the  expression  used  In  Region  II.  It  Is  accurate 
to  about  two  parts  In  10^  for  any  (x,y)  In  the  region. 

For  Region  III  a two  point  Gauss-Laguerre  quadradure  formula  gives 


w(Z)  - IZ 


®1  ^ ®3 


2-  ' 2 
Z Z -B^ 


(4.4-10) 


- 0.5124242 


Bg  - 0.05176536 


B^  - 0.2752551 


B^  - 2.724745 


Following  the  same  procedure  as  before,  the  real  portion  of  the  function 
can  be  acertalned  as 

? 2 2 2 2 2 
, ^ 2x  y-(x  -y  -B,)y  2x  y-(x  -y  -B,)y 

RfwCZ)}  - B I + B3  ! (A.A-11) 

I J ^ (x^V-B,)  + C2xy)  (3t2.y2_g  .2  ^ (2^j2 
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r 
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1 
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This  expression  la  accurate  to  about  one  part  In  10  for  Cx^y)  in 

Region  III.  The  FORTRAN  IV  program  implementing  Eqs.  4.4-4,  4.4-9,  I 

and  4.4-11  is  listed  in  Appendix  C.  i 

! 

4.5  Comparison  of  Voigt  Algorithms 

In  order  to  select  one  of  the  previous  methods  for  evaluation 

of  the  Voigt  function  it  was  necessary  to  design  computer  programs 

to  verify  the  accuracy  of  each  method  and  establish  some  relative 

times  for  execution.  Additionally  the  ease  of  Implementation  and 

program  alse  of  each  method  had  to  be  considered. 

It  waa  assumed  that  any  time  savings  in  the  computation  of  the 

parameters  x and  y would  be  as  beneficial  to  one  method  aa  another. 

Tharafore  tha  axacutlon  time  of  a method  is  defined  to  be  the  time 

that  axplraa  from  tha  issuance  of  a call  statement  to  tha  return  of 

a value  by  tha  subroutine.  Compilation  time  of  tha  program  waa  not 

included  nor  waa  any  input  or  output.  Timing  of  tha  routines  became 

a problem  and  several  methods  of  time  determination  wara  tried.  In 

tha  and  it  was  decided  to  use  tha  execution  time  returned  by  tha 

operating  syatem  of  tha  computer.  That  should  account  for  all  the 

tiaie  that  a program  has  control  of  tha  central  processing  unit  (CPU) 

and  ia  probably  as  accurate  as  any  timing  routine. 

Accuracy  datarmlnatlons  were  somewhat  simpler.  Faddayeva  and 
12 

Terantdv  have  tabulated  values  for  v(Z)  whan  o < x < 5 and  o < y < 5 
to  six  decimal  figures.  For  regions  outside  of  that  range  hrayaon's^^ 
method  was  used  to  set  the  true  value  of  w(Z) . It  was  assumed  that 
alnca  his  work  had  been  published  it  would  be  as  accurate  as  he  claimed. 
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All  programs  were  processed  by  an  IBM  360/65  computer. 


A. 5.1  Accuracy  Tests 

Three  separate  tests  were  conducted  to  determine  the  accuracy 
of  each  method.  As  an  initial  screening  Klelkopf's  method  and  the 
method  proposed  here  were  pitted  against  Drayson's  method.  A synop- 
sis of  the  results  are  shown  In  Table  A. 5-1. 

This  first  test  was  conducted  using  double  precision  operations 
to  increase  accuracy  In  all  measurements.  The  table  does  not  show 
all  the  digits  used,  which  Is  why  entries  such  as  the  one  for  x ■ 0.5, 
y " 4.0  appear  to  have  no  error  when  actually  there  Is  a slight  dif- 
ference. As  the  table  Indicates,  Klelkopf's  method  has  accuracy 

errors  that  are  relatively  large.  Ex^uIllnatlon  of  his  values  with 

12 

true,  tabulated  values  [Faddeyeva  and  Terentdv  ] revealed  unaccept- 
able margins  of  error.  Thus  the  same  conclusion  drawn  by  Drayson  In 
his  search  for  an  accurate  Voigt  Algorithm  was  reached.  Klelkopf's 
method  Is  not  accurate  enough  for  atmospheric  work  [Drayson^^] . 

The  second  test  was  designed  to  verify  the  accuracy  of  both  the 
proposed  method  and  Drayson's  method.  The  regions  used  by  the  proposed 
method  were  selected,  one  at  a time,  and  one  hundred  values  of  the 

function  were  computed  by  each  method.  The  computed  values  were 

12 

compared  to  those  tabulated  by  Faddeyeva  and  Terentdv  for  o < x < 5, 

0 < y < 5.  In  Region  III,  where  x > 5,  y > 5 the  values  returned  by 
the  two  methods  were  compared  with  each  other,  since  no  tabulated 
values  were  available. 

Figure  4.5-1  depicts  the  regions  used  by  the  proposed  method. 
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True  Value 
Olrayson's) 

Klelkopf ' s 
Value 

Proposed 

Value 

Percent  leletlve 

Error 

Klelkopf  1 Proposed 

.426050 

.426182 

.426042  I 

3.09*10'^ 

1.84* 10"^ 

.413984 

.415461 

.413989 

3.56*10"^ 

1.21 ‘lO"^ 

.391218 

.394833 

. 391236 

9. 24- 10"^ 

4. 71* 10"^ 

.254974 

.254805 

.254970 

6. 62* 10"^ 

1.56 *10"* 

.251677 

.244411 

.251692 

2. 88. 10"^ 

5. 96* 10"* 

.245279 

.224884 

.245226 

8.31* 10"^ 

2.16*10"* 

.178842 

.180084 

.178841 

6. 94- 10"^ 

4. 11* 10"* 

.177581 

.171631 

.177581 

3.35*10"^ 

2.78*10"* 

.175105 

.155897 

.175105 

1.09* 10"^ 

2.56*10"* 

.136923 

.139572 

.136924 

1.93*10"^ 

1.28*10"* 

.136329 

.132635 

.136330 

2. 70* 10"^ 

1.05*10"* 

.135155 

.119779 

.135155 

1.13* 10"^ 

6.60*10"* 

.110663 

.114381 

.110664 

3. 35 •10"^ 

2.78*10"* 

.110341 

.108533 

.110342 

1.63'10"^ 

2.20*10"* 

.109702 

.092719 

.109703 

1.55-10"^ 

1.75*10"* 

Table  4.5-1.  Initial  accuracy  tests.  Drayson's  values  used  as 
true.  (Not  all  digits  shown  In  value  columns.) 


i i 
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Tha  shaded  areas  are  the  range  of  the  points  used  for  comparison. 

The  points  were  selected  to  Include  as  much  of  the  area  of  each  region 
as  possible.  In  Region  I 

y ■ 0 to  1.8  In  Increments  of  0.2, 

X ■■  0 to  2.7  In  Increments  of  0.3. 

In  Region  II 

y - 2 

y ■ 2.5  to  4.6  In  Increments  of  0.3, 

X > 0 to  4.5  In  Increments  of  0.5, 
and  y ■■  0 to  3.5  In  Increments  of  0.5, 

X ■■  3 to  4.6  In  Increments  of  0.4. 

In  Region  III 

y - 5 to  17  In  Increments  of  3.0, 

X ■ 2 to  11  In  Increments  of  1.0, 
and  y « 0.5  to  5 In  Increments  of  0.5, 

X 5 to  21  In  Increments  of  4.0. 

The  results  of  this  test  are  summarized  In  Table  4.5-2.  The  peak 
deviation  In  Region  I for  the  proposed  method  Is  larger  than  the  peak 
for  Drayson,  but  the  actual  error  for  that  single  point  was  only  three 

4 

parts  In  10  . The  root  mean  square  error  was  well  within  acceptable 
limits. 

For  Region  II  the  peak  deviations  were  nearly  the  same,  while 
the  root  mean  square  error  for  the  proposed  method  was  better  than  Dray- 
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Region 

Drayson ' < 
Peak  Deviation 

3 Method 

RMS  Deviation 

Proposed  Method 

Peak  Deviation  RMS  Deviation 

I 

1.72504.10"^ 

4.29676*10"^ 

3.89991*10"^ 

4.07495*10"® 

11 

1.41329*10"^ 

1.41402*10"^ 

1.41342*10"® 

1.42011*10"® 

III 

*3.41969*10"^ 

*6.03836*10"® 

*3.41969*10"® 

*6.03836*10"® 

^Deviation  of  one  method  from  another. 


Table  A. 5-2.  Results  of  accuracy  test  two. 
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aon'a.  In  Region  III  the  deviation  of  one  method  from  the  other  la 
amall,  reflecting  the  two  different  q^uadrature  techniques  used. 

The  final  accuracy  test  uaed  the  regions  suggested  by  Drayson, 
one  at  a time.  Figure  4.5-2  depicts  the  coverage  of  each  area.  Rather 
than  using  tabulated  values,  the  results  obtained  with  Drayson' a 
method  were  compared  to  those  obtained  using  the  proposed  method. 

Table  4.5-3  summarises  the  results  of  the  test. 

It  Is  clearly  seen  from  Table  4.5-2  and  Table  4.5-3  that  the 
proposed  method  yields  values  of  the  Voigt  profile  that  are  accurate 
enough  for  atmospheric  work.  Root-mean-square  errors  lu  Icate  accura- 
cies that  normally  exceed  tabulated  values  of  the  profile. 


■1  ! 
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Figure  4.5-2. 


Coverage  of  Dra> 


8 10  12  14  16 

X 


Region 

Peak  Deviation 

1 

RMS  Deviation 

1 

Nundier  Points  Used 

I 

5. 27915* 10"^ 

2.46145*10"^ 

452 

11 

3. 44067* 10"^ 

1.34552*10"^ 

283 

Ilia 

6.92595*10"^ 

2.06357*10"® 

744 

Illb 

3.41969*10"^ 

1.99151*10"® 

420 

' I'l 
A 

' ^ 

Table  A. 5-3.  Deviations  of  Proposed  method  from  Drayson's  I 

method  using  Drayson's  regions.  | 
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4.5.2  Execution  Time  Testa 


Four  testa  Here  conducted  utilizing  Drayson's  method  and  the 
proposed  method  In  an  effort  to  determine  vhlch  algorithm  would 
save  the  most  computer  execution  time.  For  each  of  the  tests  a 
small  main  routine  was  written  that  fixed  the  values  of  x and  y and 
then  called  the  routine  being  tested.  The  master  program  of  some 
tests  performed  Input/output  or  accuracy  measurements.  The  other 
main  routines  performed  no  extraneous  operations.  The  tests  Involved 
merely  replacing  that  portion  of  the  program  Implementing  Drayson's 
method  with  the  proposed  method  emd  resubmitting  the  job  for  execution. 

The  first  test  used  the  regions  of  the  proposed  method,  selecting 
points  from  one  region  at  a time.  This  Is  a slightly  biased  test 
since  the  proposed  algorithm  remained  In  one  region  while  Drayson's 
algorithm  jumped  from  region  to  region,  depending  on  the  value  of 
the  parameters.  However  Drayson  contended  that  his  method  was  the 
most  rapid  available,  and  placed  no  stipulations  on  how  parameters 
were  selected.  The  main  routine  for  this  test  performed  Input  and 
output  as  well  as  accuracy  measurements  for  both  methods.  Therefore, 
the  times  were  somewhat  over  true  execution  time,  but  reflected  the 
relative  speed  of  each  method.  Table  4.5-4  lists  the  results  of  this 
test.  The  method  proposed  In  this  work  reduced  the  time  to  approximate 
the  Voigt  profile  by  an  average  of  26. 6Z. 

A second  test  was  conducted  using  the  regions  proposed  by  Dray- 
son, selecting  points  from  one  region  at  a time.  This  has  the  same 
bias  as  the  first  test  except  that  It  Is  In  Drayson's  favor  this  time. 


Region 

Drayson's  Method 

CPU  Seconds  % 

Proposed  Method 

CPU  Seconds  X 

I 

4.05 

100 

2.92 

72.09 

II 

4.14 

100 

3.00 

72.46 

III 

3.87 

100 

2.92 

75.45 

Table  4.5-4.  Results  of  time  testa  conducted  using  regions 
and  points  indicated  in  Fig.  4.5-1. 
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One  hundred  points  were  selected  from  each  region.  The  results  of 
the  test  are  shown  In  Table  4.5-5.  The  main  routine  for  this  test 
perfoirmed  output  operations  as  well  as  calling  the  method  being  tested. 
The  mechod  proposed  by  this  work  reduced  the  time  to  approximate  the 
Voigt  profile  by  an  average  of  27.5%. 

An  offshoot  of  the  accuracy  test  conducted  using  the  regions  and 
points  shown  In  Fig.  4.5-2  was  that  the  times  for  execution  were  avail- 
able. This  comprised  the  third  speed  test.  The  main  routine  In  this 
case  performed  output  and  deviation  analysis  in  addition  to  calling 
tha  Voigt  algorithm.  Again  the  times  were  somewhat  In  excess  of  true 
execution  time,  but  still  reflected  the  relative  speed  of  each  method. 
Table  4.5-6  summarizes  this  test.  The  proposed  method  reduced  the 
time  to  evaluate  the  Voigt  profile  by  an  average  of  29.5%. 

A final  test  was  conducted  that  did  not  use  regions.  Ten-thousand 
points  were  selected  as  follows 

y > 0 to  19.8  In  Increments  of  0.2, 

X ••  0 to  24.75  In  Increments  of  0.25. 

Tha  master  routine  performed  no  Input  or  output  and  no  other  calcula- 
tions. This  was  strictly  comparing  the  time  for  Drayson's  method  and 
the  proposed  method  to  evaluate  ten-thousand  points.  The  results  are 
enumerated  In  Table  4.5-7. 

4.5.3  The  Selected  Voigt  Algorithm 

In  view  of  the  results  of  the  time  tests  conducted  the  selection 
was  not  difficult.  Both  methods  yield  the  accuracy  necessary  for  good 
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Region 

Dra^son's  Method 
CPU  Seconds  Z 

Proposed  Method 

CPU  Seconds  X 

I 

3.30 

100 

2.37 

71.81 

II 

3.22 

100 

2.52 

78.26 

Ilia 

3.90 

100 

2.52 

64.61 

mb 

3.59 

100 

2.70 

75.21 

Tabl*  4.5-*5.  Results  of  time  tests  conducted  using  100  points 
from  each  region  Indicated  In  Fig.  4.2-1. 
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Region 

Drsyson's  Method 
CPU  Seconds  X 

Proposed  Method 

CPU  Seconds  X 

I 

3.45 

100 

2.29 

66.38 

II 

3.12 

100 

2.49 

79.80 

Ilia 

3.72 

100 

2.69 

72.31 

Illb 

3.69 

100 

2.42 

65.58 

Y 


Tabl«  4.5-6.  Results  of  time  tests  conducted  using  regions  snd 
points  Indicated  In  Fig.  4.5-2. 

t’ 
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atmo8ph*rlc  calculation  wd  both  methods  are  easy  to  program.  Drayson's 
algorithm  contains  some  64  source  statements  while  the  proposed  method 
has  40.  Drayson's  method  was  run  In  single  precision  and  the  proposed 
method  In  double  precision.  It  is  felt  that  the  proposed  algorithm 
could  be  executed  In  single  precision  and  still  retain  the  accuracy 
required.  The  conversion  to  single  precision  would  only  increase  the 
speed  of  the  method.  Therefore  it  is  the  proposed  method  that  will 
be  used  to  evaluate  the  Voigt  profile  In  the  transmittance  differential 
equation. 


CHAPTER  V 

ABSORPTION  LINE  PARAMETERS  AND  UNITS 

5.1  Preface 

This  chapter  brings  together  the  remaining  topics  that  are 
necessary  to  develop  the  computer  solution  to  the  transmittance  dif- 
ferential equation.  A brief  review  of  the  equations  to  be  used  Is 
given  and  the  known  and  unknown  quantities  appearing  in  the  relations 
are  Identified.  The  Input  data  necessary  to  solve  the  equations  are 
also  identified.  The  value  of  physical  constants  used  in  the  calcu- 
lations is  given  at  this  time.  The  mixing  ratio  M(P)  Introduced  In 
Eq.  3.2-11  is  derived  and  the  combined  extinction  coefficient 
relation  Is  developed.  Then  the  units  of  all  the  relations  used  In 
this  work  are  dimensionally  analyzed  to  assure  they  are  correct. 
Additionally  some  unit  conversion  factors  are  Introduced.  The  data 
available  as  atmospheric  absorption  line  parameters  are  presented 
together  with  some  general  remarks  on  their  derivation  and  compilation. 
Finally  a computer  routine  that  is  used  to  handle  the  tremendous 
amount  of  data  available  Is  developed  and  explained. 

5 . 2 Review  of  Transmittance  Equations 


The  absorption  coefficient  relations 


dt 

2 2 
+ (x-t)'^ 


(5.2-2) 


(5.2-3) 


and  the  line  strength  function 


'ij 


(5.2-4) 


For  the  purpose  of  dimensional  analysis  the  summations  over  1 and  j 
are  dropped.  The  variables  appearing  in  the  above  equations  must 
also  be  defined, 

y = — (ln2)^^^  . (5.2-5) 

“d 


v-v  , 

5.  (in2)^'^  , (5.2-6) 

“d 


(5.2-7) 
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It  is  seen  that  in  order  to  solve  these  equations  it  is  necessary  to 
specify  : 

S^...  a reference  line  strength, 

Q ...  the  vibrational  partition  function, 

T^...  the  reference  temperature, 

T ...  the  temperature  at  which  transmittance  is  to  be  calculated, 
c ...  the  exponent  that  expresses  the  rotational  partition  function, 
E ...  the  energy  of  the  line, 
k ...  Boltzman's  constant, 

M ...  The  molecular  weight  of  the  absorber, 
v^...  the  center  wavenumber  of  the  line, 
a^...  the  reference  Lorentzlan  line  half-width, 

P^...  the  reference  pressure, 

F ...  the  pressure  at  which  transmittance  is  to  be  calculated, 

V ...  the  wavenumber  at  which  transmittance  is  to  be  calculated, 

M(P) . the  mixing  ratio  of  the  absorber. 

These  variables  may  be  grouped  into  four  broad  catagorles.  Those 
dependent  upon  the  absorbing  line  (S^,  E,  v^,  » those  dependent 

upon  the  absorber  type  or  gas  (Q,  c,M,M(P)),  those  dependent  upon 
the  atmosphere  (T,  P) , and  finally  those  that  are  constant  (T^,  P^, 
k,  v) . Actually  v is  not  constant  since  if  the  transmittance  at 
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another  wavenmber  is  desired,  the  value  of  v must  be  changed.  Data 
must  be  obtained  then,  from  the  atmosphere,  from  the  abosrblng  gas, 
and  from  the  absorbing  line. 

5.3  The  Mixing  Ratio  and  The  Combined  Extinction  Coefficient 

The  mixing  ratio  M(^)  Introduced  In  chapter  three  expresses  the 
relation  between  total  pressure  and  the  pressure  of  only  one  gas 
Cabsorber)  under  consideration.  In  the  development  of  the  Lambert 
Beer  law  (Eq.  3.2-7)  it  was  stated  that  the  reduction  in  Intensity 
of  radiation  was  proportional  to  the  amount  of  gas  present  In  the 
path.  If  u represents  the  gas  amount,  then  du  Is  the  Incremental 
gas  amount  and  the  Lamber-Beer  Law  may  be  written 

T = exp (-Kp (x) dx)  = exp(-Kdu).  (5.3-1) 

So  that  the  gas  amount  Is 

du  = p (x)  dx.  (gm/cm^)  (5.3-2) 

gas 

If  the  right  side  of  Eq.  5.3-2  is  divided  by  the  density  of  the  gas 
at  standard  temperature  and  pressure  (STP) , the  units  of  the  gas 
amount  are  changed  to  centimeter-atmospheres  (cm-atm) . This  unit  is 
more  commonly  used  In  atmospheric  work  since  It  Is  easier  to  compute 
the  equivalent  amount  of  an  absorber  In  a path  length  as 

AL  = cRP  (5.3-3) 

where  c Is  the  fractional  concentration  of  the  absorber,  R the  path 
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I 


i 

[ 

1 

[ 

f 


length  (cm),  and  P the  pressure  (atm).  Converting  Eq.  5.3-2  to 
cm-atm  yields 


du 


P 

- - dx 

^gas(STP) 


(cm-atm) . 


In  order  to  Introduce  a mixing  ratio  multiply  and  divide  the  relation 
by  the  density  of  air. 


du  = — — dx  (cm-atm) 

•^alr  ^gas(STP) 

^alr  ^eas 

du  = dx  (cm-atm) 

^gas(STP)  ^alr 


Now  define  the  density  mixing  ratio 


M,  = ^ ^ . (5.3.4) 

*^air  gm/cm 

where  the  densities  are  at  some  temperature  and  pressure.  Now  gas 
amount  may  be  written 


du  = M,  dx.  (cm-atm)  (5.3-5) 

^gas(STP)  ^ 

The  relation  between  the  volumetric  mixing  ratio  measured  In  parts 

per  million  (ppm)  and  the  density  mixing  ratio  (M^)  is  given  by 
49 

[Weast  ] as 
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M • 10  • Mt 


(5.3-6) 


where  Mt  Is  the  molecular  weight  of  air  and  gas  respectively.  Now 


p ^ M • Mt  *10 
air  ppm  gas 


‘’gas(STP)  “’^alr 


M • 10  Mt  p . 
du  = dx. 

‘^gasCSTP)  air 


(5.3-7) 


The  hydrostatic  approximation  given  by  Eq.  2.2-1  relates  density  to 
pressure  In  a manner  desired. 


j dP 
p , dx  = — 
air  g 


(5.3-8) 


Here  F Is  pressure  and  g.  Is  the  acceleration  due  to  gravity.  So 
finally 


M • 10  Mt  dP 

du  = —222 — \ . E.^1. (cm-atm) 

^gas(STP)  “^alr  ® 


(5.3-9) 


du  - 


M^  dP 
d 

^gas(STP)  g 


(cm-atm) , 


(5.3-10) 


Either  Eq.  5.3—9  or  5.3—10  may  be  used  to  calculate  the  gas  amount. 
The  selection  of  one  over  the  other  Is  dependent  upon  the  gas  and  the 
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mixing  ratio.  For  most  atmospheric  gases  Eq.  5.3-9  Is  used.  In  the 
case  of  water  vapor  Eq.  5.3-10  Is  used.  The  calculated  gas  amount 
la  multiplied  by  the  extinction  coefficient  In  the  calculation  of 
transmittance.  As  might  be  expected  there  Is  a different  extinction 
coefficient  (K)  and  gas  amount  (du)  for  each  absorber  (gas)  In  the 
atmosphere. 

When  more  than  one  absorber  (gas)  Is  to  be  considered  In 
atmospheric  transmittance  the  Index  of  the  exponential  must  be  the 
sum  of  all  optical  depths. 

T exp(-|  KjMjdP). 

Since  every  absorber  has  a different  mixing  ratio  and  extinction 
coefficient,  each  must  be  computed  and  multiplied  before  they  can 
be  added. 

Table  5.3-1  lists  the  mixing  ratios  In  parts  per  million  (ppm) 
for  those  atmospheric  absorbers  whose  fractional  concentrations  do 
not  vary  as  a function  of  altitude.  It  also  lists  the  molecular 
weights  of  the  gases  of  constant  mixing  ratio.  Table  5.3-2  gives 
the  density  at  STP  of  the  atmospheric  gases  that  absorb  Infrared 
radiation. 

For  cases  where  the  mixing  ratio  Is  not  constant  (water  and 
ozone)  a different  value  of  M(P)  must  be  used  at  each  altitude 
where  Transmittance  Is  calculated.  Table  5.3-3  lists  the  mixing 
ratios  for  the  variable  atmospheric  constituents  water  and  ozone  at 
different  altitudes. 

For  the  atmosphere  only  seven  absorbers  need  be  considered. 
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Atmospheric 

Constituent 

Chemical 

Formula 

Parts  per  Million 
by  Volume 

Molecular 

Weight 

Air 

10' 

28.97 

Carbon  Dioxide 

CO2 

330 

44 

Nitrous  Oxide 

N2O 

0.28 

44 

Carbon  Monoxide 

CO 

0.075 

28 

Methane 

CH^ 

1.6 

16 

Oxygen 

°2 

2.095-10^ 

32 

Table  5.3-1.  Mixing  Ratios  and  Molecular  Weights  of  Non-Varying 

31 

Atmospheric  Constituents.  [McClatchey,  et  al.  ] 
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Atmospheric 

Constituent 

Chemical 

Formula 

Density  at  STP 
(gm/ca^) 

Carbon  Dioxide 

CO2 

1.9768441* 10"^ 

Nitrous  Oxide 

N2O 

1.9781371*10"^ 

Carbon  Monoxide 

CO 

1.2510101*10"^ 

Methane 

CH4 

7.1678376*10"^ 

Oxygen 

°2 

i 

1.4286545*10"^ 

Ozone 

O3 

2.1436282*10"^ 

Water 

H2O 

1.0000010 

Table  5.3-2  Density  of  Atmospheric  Absorbers  at  Standard  Temperature 


36 

and  Pressure.  [Nelson  ] 


Height 

(Rm) 

Temperature 

("K) 

Mixing  Ratio  (M^) 

Water  Vapor  Ozone 

0 

294 

1.1754* 10"^ 

5.0378-10"® 

5 

267 

1.3867- 10"^ 

9.1527-10"* 

10 

23'^ 

1.5388- 10"^ 

2.1639*10"^ 

15 

216 

3.6122* 10"^ 

9.0304*10"^ 

20 

218 

4.7604-10"® 

3.5967-10"® 

25 

224 

1.5625-10"^ 

6.9962-10"® 

30 

234 

2.7231-10"^ 

1.5128-10"^ 

35 

245 

1.6873-10"^ 

1.4112-10"^ 

40 

258 

1.2913-10"^ 

1.2312*10"^ 

45 

270 

1.0814-10"^ 

7.3989-10"® 

50 

276 

6.6232-10"® 

4.5204-10"® 

70 

218 

2.0876-10"® 

1.2824-10"® 

100 

210 

2.0-10"® 

8.6-10"® 

Table  5.3-3.  Mixing  Ratios  for  Water  Vapor  and  Ozone  from  A Midlatitude 
Summer  Atmosphere.  [McClatchey,  et  al.^^1 
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Five  of  these  have  a constant  mixing  ratio  and  two  have  variable 
mixing  ratios.  In  any  case  we  are  concerned  with  the  product  K du. 
Define  the  combined  extinction  coefficient 


ABSK  = Z K du 

j J J 


(5.3-11) 


Let  the  value  of  the  subscript  j vary  from  one  to  seven  to  account  for 
all  gases.  Further,  establish  the  following  relations 


1 . 
2 . 

3 . 

4 . 

5 . 

6 . 
7 ., 


. denotes  water  vapor, 

. denotes  carbon  dioxide, 

. denotes  ozone, 

. denotes  nitrous  oxide, 

. denotes  carbon  monoxide, 
. denotes  methane, 

. denotes  oxygen. 


(5.3-12) 


Kj^dUj^  is  the  extinction  coefficient  for  water  vapor,  K2du2  for  carbon 
dioxide,  and  so  on.  Defining  the  variable  L as  the  counter  for  the 
level  of  the  atmosphere  we  can  write 


ABSK(L)  = Kj(L)dUj(L) 


The  gas  amount  du^ (L)  will  have  seven  different  expressions;  one  for 


each  gas.  For  water  vapor  and  ozone  the  gas  amount  will  be  a function 
of  atmospheric  level  and  can  be  written 


du(L)  = — 


Md(L) 


g P 


gas(STP) 


(5.3-13) 


For  the  other  gases  the  gas  amount  is  not  so  dependent  upon  atmospheric 
level.  It  is  only  a function  of  pressure  (which  is  dependent  upon 
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level) 


M 


duCL) 


PP”> 


10"^  Mt 


gas 


dP. 


^gas(STP) 


(5.3-14) 


Thus  the  combined  extinction  coefficient  can  be  written  as 


ABSK(L)  = 


+ K (L)  (330  ■ 10-^)  (44) 

(1.97  • 10"^) (28.97) 


Kott)  , + V (L)  _,X-28  • .. 

« - * - -I  J ^ fH  /X-t  - 


2.14  • 10 


(1.97  • 10  ■')  (28.97) 


K5(L)  ~ • 10'^  (28)  ^ (1.6  . 10-^(16)  ^ 

(1.25  . 10“^) (28.97)  ^ (7.16  • lO'S (28.97) 


K7(L) 


(2.095  • 10~^)(32)  I 
(1.42- . 10"^)(28.97)J 


where  the  value  of  Mpp^.Mtg^^.Mt^^^and  ^ Inserted 

2 

into  the  equations.  Assuming  a value  of  980.665  cm/sec  for  g yields 
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ABSK(L)  - [1.02  • 10"^K^(L)M^(L)  +2.59  * 10"'*K2(L)  + 

4.76  • 10“^K-(L)M  (L)  + 2.19*10"V(L)  + 

^ ^ ^ (5.3-15) 

5.9  * lO’^Kjd)  +1.25  • lO'^K^a)  + 

1.65  • 10~^K^(L)]dp 

This  relation  Is  used  to  obtain  the  combined  extinction  coefficient 
at  each  level  (L)  of  the  atmosphere. 

5.4  Units 

In  the  calculation  of  transmittance  It  Is  necessary  to  determine 
gas  amounts  and  extinction  coefficients.  For  different  applications 
or  when  different  data  Is  available.  It  may  be  necessary  to  convert 
units.  Therefore  a brief  treatise  Is  given  on  conversion  factors. 

According  to  Avogadro's  hypothesis  the  molecular  weight  (M) , In 
gr2uns,  of  any  gas  occupies  22.4  liters  at  standard  temperature  and 
pressure  (STF) . Hence  one  cubic  centimeter  of  gas  at  STP  weighs 
M 

T—  grams.  Since  one  centimeter-atmosphere  at  STP  Is  equlva- 

2.24-10  2 

lent  to  a length  of  one  centimeter  of  gas  at  STP  per  cm^.  It  Is 

possible  to  write  the  following  relation 


1 cm-atm  of  gas 


M 

2.24-10^ 


2 

gm-cm  , 


(5.4-1) 


where  M Is  molecular  weight. 

3 1 

Also  since  one  cm  of  gas  at  STP  contains  t 

2.24-10 

one  mole  of  gas  contains  Avogadro's  number  of  molecules 

3 

possible  to  write  that  one  cm  of  gas  at  STP  contains  — 


moles  and 


, It  Is 

23 

6.02-10'^^ 

2.24-10^ 
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molecules , 


l(cm-atm)g^p  ^ 2.69*10^^  molecules/cm^ 


(5.4-2) 


Equation  5.4-1  is  valid  Independent  of  the  gas.  Since  for  water 
vapor  it  is  possible  to  write  1 gm-cm  in  terms  of  (cm-atm)g^p  it 
is  shown  for  water  vapor  that 

1 precipitable  centimeter  (H20)=3. 34* 10^^  (molecules/cm^)  (5.4-3) 

2 

Thus  the  unit  molecule/cm  is  independent  of  the  nature  of  the 
absorbing  gas  and  basic  to  all  gases. 


Some  other  relations  that  have  proven  useful  in  this  work  are; 

. -Iv  frequency  (1/sec)  , 

wavenumber  (cm  ) - speed  of  lighV  Tcm/secT ' 


(5.4-4) 


energy (joules)  = energy  (wavenumbers) 'Planck's  constant 
(joules-sec) "Speed  of  light  (cm/sec) 


E(j)  = [E(cm  ■'■)  ] [h(j-sec)  ] [c(cm/8ec)  ] 


(5.4-5) 


5 2 

1 atm  ® 1013.25  milibars  = 1.01325  ' 10  newtons/m  (5.4-6) 

Table  5.4-1  lists  the  physical  constants  that  were  used  in  this  work. 


5.4.1  Dimensional  Analysis  of  Transmittance  Equations 

From  the  equations  used  in  this  work  it  can  be  seen  that 


must  be  dimensionless. 


Equation  5.2-2  gives 


+ (x-t)^ 


(5.4-7) 
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Physical  Constant 


Symbol 


Valus 


Unit 


Acceleration  of  gravity 

g 

980.665 

cm/ssc^ 

Planck  constant 

h 

6. 626196 *10“^^ 

Joule-sec 

Speed  of  light 

c 

2.997925*10^° 

cm/sec 

Avogadro  constant 

Me 

6.022169*10^^ 

molecule  ^ 

doltsman  constant 

k 

1.380622*10'^^ 

Joule/*K 

Reference  Pressure 

P 

o 

1013.25 

mlllbars 

Reference  Temperature 

T 

o 

296 

•k 

PI 

Tt 

3.141592654 

Table  5.4-1.  Physical  constants  used  In  this  work.  (Adapted  from  Weast  } 


I- 


with 


V 


— and  x 

“d 


v-v 

o 


(ln2) 


1/2 


Provided  that  ,a  ,v,v  ,1  , c “1  4 ^ 1 

L T)  o all  have  units  of  cm  the  integral 

portion  of  Eq.  5.4-7  is  dimensionless.  Therefore  K will  have  units 

of  given  by 


K 

o 


in.: 


a/2 


(5.4-8) 


where 


(5.4-9) 


As  long  as  and  T have  the  same  units  their  ratio  is  dimensionless. 
The  partition  function  Q has  a numerical  value  and  is  dimensionless 
and  the  index  of  the  exponential  is  dimensionally 

(Joules) f ^ ^ \ 

(Joules/'K)  ('’K)  I °K  | » 

provided  the  energy  E is  given  in  joules.  If  it  is  not  it  must  be 

converted.  (See  Eq.  5.4-5).  Therefore  the  strength  S will  have 

units  of  S . 

o 


r.  c d 

S = S = 

o 


cm 


molecule  - cm 


-2 


(5.4-10) 
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From  Eq.  5.4-8  it  can  be  seen  that  K , and  therefore  K will  have 

o 

units  of 


cm 


-1 


molecule  - cm 


-2 


cm 


-1 


2 

cm 


molecules 


(5.4-11) 


From  the  discussion  of  section  5.3  recall  that  the  quantity 


of  Eq.  5.4-7  Is  really  gas  amount  du  and  Is  given  by  either 

g 

Eq.  5.3-9  or  5.3-10.  In  either  case  It  has  units  of  cm-atm.  There- 
fore 

7 

• cm-atm 


KMdP  d cm' 
g 


molecule 


Utilizing  Eq.  5.4-2  It  Is  possible  to  make  the  above  relation 
dimensionless . 

19 

, 2.69  • 10  molecules 

1 3 

cm  - atm 


KMdP  d 
g 


cm 


molecules 


• cm-atm 


19 

2.69  • 10  molecules 


cm  - atm 


Thus,  all  calculations  will  be  dimensionally  correct  provided 

(1)  Ojj,v,v^  are  in  wavenumbers  (cm  ^) , 

_-l 

(2)  S has  units  of 


cm 


molecule-cm 

(3)  E has  units  of  Joules, 


-2 


; i 
; 1 
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(4)  The  extinction  coefficient  K Is  multiplied  by  2.69*10^®  molecules ^ 

cm  -atm 


C5)  Eq.  5.2-1  is  used  to  calculate  transmittance  where 


M 


10"^  • Mt 


gas 


M (P)  . _E2“ 

J Pgas(STP)  «^lr 


or 


M. 


M (p)  ii 

’^gasCSTP) 


(6)  The  physical  constants  of  Table  5.4-1  are  used  where  needed. 


5.5  The  Absorption  Line  Parameters 

About  10  years  ago  various  Individuals  began  to  compile  spectro- 
scopic data  on  the  vibrational-rotational  lines  of  atmospheric  gases. 

This  continued  mostly  as  Individual  work  until  about  5 years  ago  trfien 
an  effort  at  the  Air  Force  Cambridge  Research  Laboratories  (AFCBL)  was 
made  to  continue  the  data  collection  with  the  aim  of  providing  a complete 
set  of  data  for  all  vibrational-rotational  lines  of  naturally  securing 
molecules  of  significance  In  the  terrestrial  atmosphere.  With  such 
data  at  hand  It  would  be  possible  to  compute  the  transmittance  appro- 
priate for  atmospheric  paths.  The  work  at  AFCRL  Includes  data  on  water 
vapor,  carbon  dioxide,  osone,  nitrous  oxide,  carbon  monoxide,  methane, 
and  oxygen  to  date.  All  of  these  molecules  except  oxygen  are  minor 
constituents  of  the  atmosphere,  but  nonetheless  represent  most  of  the 
abaorptlon  lines  in  the  visible  and  Infrared. 

In  order  to  compute  transmittance  due  to  a given  spectral  line 


1 
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in  th*  atsoaphere  It  la  nacaaaary  to  daacrlba  tha  ahaorptlon  coafflciant 
aa  a function  of  fraquency  of  aacK  llna.  Tha  four  aaaantial  parnatara 
for  aach  llna  ara  tha  raaonanca  fraquancy  v^,  tha  Intanalty  par  abaorh- 
Ing  nolacula  S,  tha  Lorant*  llna  width  paranatar  o^,  and  tha  anargy  of 
tha  lowar  atata  B.  Tha  fraquancy,  v^,  la  Indapandant  of  both  tampara* 
tura  and  praaaura.  Tha  intanalty,  S,  la  praaaura  Indpandant  and  Ita 
taaparatura  dapandanca  can  ba  calculatad  from  v and  E. 

PTlvatlon  of  Llna  Paramatars 

Tha  four  paramatara  E,  S,  and  muat  of  couraa  ba  darlvad 

from  axparlmantal  obaarvatlona,  aubjactad  to  data  raductlon  In  tha 
framawork  of  tha  ganaral  thaorlaa  of  molacular  apactroacopy.  Tha 
complaxlty  naadad  to  approach  tha  problam  dapanda  both  on  tha  typa  of 
iwlacula  and  tha  accuracy  of  obaarvatlonal  data.  In  tha  darlvatlon  of 
tha  paramatara,  aquations  and  mathoda  vara  u-ad  for  llnaar  trlatomlc 
and  diatomic  molaculaa  (CO2,  N^O,  CO)  and  nonllnaar  trlatomlc  molaculaa 
(HjO,  Oj).  Mathana,  CH^,  a spharlcal  top,  la  a spaclal  caaa,  aa  la 
tha  diatomic  O^.  Tha  anargy  statas  and  tha  transition  probabllltlas 
batwaan  anargy  statas  of  tha  molaculaa  ara  daflnad  primarily  by  thalr 
numarlcal  valuas,  as  astabllshad  by  azparlmants  and  tha  quantum  numbara 
which  Idantlfy  tham.  In  naarly  all  tha  caaas  of  Intarast  mathamatlcal 
ralatlons  of  graatar  or  lassar  complaxlty  ralata  tha  numarlcal  propar- 
tlaa  to  tha  quantum  numbara. 

For  both  llnaar  and  nonllnaar  molaculaa  tha  vibrational  and  rota- 
tional statas  may  ba  charactsrlzad  by  quantum  ntaabars.  To  calculata 
tha  puraly  vibrational  part  of  tha  anargy  of  llnaar  molaculaa  it  is 
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necessary  first  to  compute  an  unperturbed  energy,  A similar 

procedure  la  used  for  nonlinear  molecules  to  compute  Then  the 

effects  of  reson2mce  perturbations  are  Incorporated  to  obtain  the  mole- 
cular vibrational  energies,  G^.  This  method  was  used  to  generate  those 
energy  levels  which  have  not  been  observed;  for  all  observed  states 
the  experimental  values  were  used. 

For  linear  molecules  the  rotational  energy  of  each  vibrational 
state  Is  a function  of  the  quantum  numbers  J and  1,  and  the  vibrational 
energy,  G^.  It  Is  given  by 

- G^+B^lJ(J+l)-£^]-D^tJ(J+l)-f^]^+H^tJ(J+l)-f^]^+....  (5.5-1) 
where  B , D , H are  constants  for  each  vibrational  state  that  are 

V V V 

either  determined  by  observation  or  calculated.  The  equations  for  B^ 
and  D are  similar  to  those  for  G and  likewise  require  modification 
through  resonance  perturbation. 

For  nonlinear  molecules  the  rotational  levels  are  calculated  from 
complicated  functions  of  the  three  different  moments  of  Inertia  of  the 
molecule  model  [Hertzberg^^] . Having  determined  the  energy  states  of 
the  different  molecules  and  their  Isotopes,  the  line  frequencies  (v^) 
are  determined  by  taking  the  differences  corresponding  to  all  allowed 
transitions.  These  depend  upon  the  selection  rules  for  the  molecule. 

Half-widths  (a^)  of  lines  are  taken  from  observed  values  or  com- 

31 

puted  for  each  molecular  species.  McClatchey,  et  al.  discusses  the 
procedure  used  to  calculate  the  half  width  of  each  gas  type.  In  cases 
where  there  Is  Insufficient  data  to  warrent  Inclusion  of  variable  half 
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1 

widths,  a constant  value  was  used. 

i 

The  intensities  of  the  lines  are  calculated  as  a function  of  fre- 
quency Cv) , the  vibrational  intensity  of  a nonrotating  molecule  (S^) , 

the  rotational  intensity  of  a rigid  nonvibrating  molecule  (SpQ,j.) , and  ! 

I 

a factor,  F,  that  takes  into  account  that  both  forms  of  motion  occur 

1 

simultaneously.  The  respresentatlons  of  F are  quite  complex  and  very 

dependent  upon  molecular  type  and  Isotope.  A complete  discussion  of 

31  ^ 

the  derivation  of  parameters  can  be  found  in  McClatchey,  et  al.  . 

In  order  to  establish  a minimum  intensity  value,  an  extreme  atmo- 
spheric path  was  considered,  assuming  gas  concentration  specified  in 
Table  5.3-1.  This  path  was  tangent  to  the  earth's  surface  and  extended 
from  space  to  space.  Lines  yielding  less  than  10  percent  absorption 
at  the  line  center  would  normally  be  omitted.  This  absolute  cutoff 
was  not  maintained  in  regions  of  strong  absorption  (relatively  weak 
lines  were  dropped)  and  regions  of  weak  absorption  (most  lines  were 
kept).  Table  5.5-1  gives  the  mean  half-widths  and  minimum  intensities 
for  each  molecule  used  in  this  work. 

Utlllalng  the  methods  outlined  above  and  the  experimental  results 

31 

of  many  other  scientists,  R.  A.  McClatchey,  et  al.  , has  complied  a 
data  tape  that  contains  the  parameters  needed  to  compute  transmittance. 

5.5.2  The  Line  Parameter  Data  Tape 

The  data  compilation  thus  far  described  resulted  in  the  accumula- 
tion of  the  four  quantities  v , S,  E,  a for  108,000  spectral  lines 
between  1 ym  and  the  far  infrared  for  the  molecular  species  listed  in 
Table  5.5-1.  Additional  identifying  information  has  also  been  supplied 
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Molecule 

Identifier 

Intensity  Minimum 

Mean  Half-Width 

HjO 

1 

3 10"^^ 

not  given 

COj 

2 

2.2  10"^^ 

0.07 

°3 

3 

3.5  10"^^ 

0.11 

N2O 

4 

3.0  10'^^ 

0.08 

CO 

5 

8.3  10"^^ 

0.06 

CH^ 

6 

3.3  10"^^ 

0.055 

°2 

7 

3.7  10"^° 

0.06 

Table  5.5-1. 


Mean  ha^f-vldths  and  Intensity  mlnlmums  used  by  McClatcbey, 
et  al.  In  compilation  of  AFCRL  data  tape. 


1 

for  each,  line  aa  indicated  In  Table  5.5-2.  A standard  con^uter  format 

was  adopted  for  card  or  card-image  Input.  The  units  of  the  data  are 

given  In  Table  5.5-2.  The  rotational  and  vibrational  Identifiers,  the 

31 

date,  and  the  Isotope  are  explained  by  McClatchey,  et  al.  , but  not 
used  here.  The  molecule  (columns  78-80)  correspond  to  the  Identifier 
of  Table  5.5-1. 

The  data  are  frequency  ordered  on  magnetic  tape  and  are  contained 
In  records  each  of  which  contain  40  card  Images  as  shown  In  Table  5.5-2. 

Each  record  Is  preceded  by  a variable  that  Indicates  the  number  of 
card  Images  within  the  record  (In  all  cases  40) . It  Is  necessary  to 
decode  each  record  from  the  tape  according  to  the  format 

no,  40(F10.3,  E10.3,  F5.3,  F10.3,  5A6,  A5,  13,  lA,  13)  23A30. 

The  last  format  (23A30)  Is  to  read  the  space  that  separates  each  re- 
cord. End  of  file  markers  are  placed  on  the  tape  at  the  following 
points:  500,  1000,  2000,  5000,  7500,  10000  cm~^.  This  results  In 
the  tape  having  seven  files. 

i 

f 5.6  The  DATSET  Computer  Subroutine 

I In  order  to  utilise  the  data  tape  of  McClatchey  It  was  necessary 

I to  write  a computer  subroutine  that  would  transfer  that  portion  of  the 

I 

data  needed  from  the  tape  to  computer  memory  or  If  desired  to  disk.  The 
subroutine  Is  used  In  conjuctlon  with  the  transmittance  program.  Some 
Job  control  language  (JCL)  Is  needed  with  the  routine  to  Identify  the 
different  files  of  the  tape  and  establish  the  disk  space  If  that  Is  to 
be  used.  The  flow  of  the  routine  Is  outlined  below  and  a listing  of 
the  program  Is  given  In  Appendix  D. 
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Data  element 

Units 

Card  Column 

Format 

V 

o 

-1 

cm 

1-10 

F10.3 

-1 

cm 

1 1 

VIA  % 

s 

0 

molecule-cm^ 

ciXU  s 

1 

a 

0 

cm”^/atm 

21.25 

F5.3 

E 

cm”^ 

26-35 

F10.3 

Rotational-Vibration  ID 

36-70 

5A6,A5 

Date 

71-73 

13 

Isotope 

74-77 

14 

Molecule 

78-80 

13 

Table  5.5-2.  Cemiputer  foraat  of  each  spectral  line’s  data.  [McClatchey, 
et  al.^^] 
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The  subroutine  has  the  following  parameters: 


WAVLO 


WAVUP 


IC!IT(T) 


MAX 


SPECS (3000,4) 


KIND (3000) 


LTOTAL 


. . . . lower  wavenumber  of  data  to  be 
extracted  from  data  tape  (lower 
limit  of  channel) 

....upper  wavenumber  of  data  to  be 
extracted  from  tape  (upper 
limit  of  channel) 

. . . .Array  containing  the  number  of 
spectral  lines  between  wavlo 
and  wavup,  by  gas  type. 

(1-  H^O,  2*  CO2,  etc.) 

. . . .Maximum  number  of  lines  that 
can  be  transfered  from  tape 
(3000) 

. . . .Array  containing  data  for  each 
line  (v^,  S^,  E,  a^) 

. . . .Array  containing  molecular 
Identifier  of  each  line 

....Total  number  of  lines  between 
wavlo  and  wavup 


. . . . Input 


. . . . Input 


.... output 


. . . .input 


.... output 


. . . .output 


. . . .output 
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There  are  error  messages  that  will  inform  the  user  If  the  program 


does  not  execute  properly.  The  routine  can  be  used  to  transfer  data 
to  disk,  (see  comment  card  In  program  listing) . The  JCL  necessary  to 
define  the  files  must  be  Included  before  using  disk. 


CHAPTER  VI 


THE  TRANSMITTANCE  PROGRAM 


6.1  Preface 

This  chapter  presents  the  object  of  the  previous  five  chapters, 
a transmittance  program.  Before  the  routine  is  presented  the  vi- 
brational partition  function  (Q)  is  developed,  the  Interatlve  solution 
used  to  solve  the  transmittance  differential  equation  is  explained, 
and  the  methods  used  to  compute  monochromatic  and  averaged  transmit- 
tance are  given.  These  final  details  are  incorporated  with  the 
equations  thus  far  developed  to  yield  the  transmittance  program.  The 
chapter  concludes  with  the  results  of  some  transmittance  calculations 
in  the  fifteen  micrometer  band. 


6.2  The  Rotational  and  Vibrational  Partition  Functions 
The  equation  for  line  strength  may  be  given  as 


S = Qr  exp 


(6.2-1) 


The  line  strength  is  dependent  upon  the  rotational  and  vibrational 

modes  of  the  molecule  in  a characteristic  manner.  The  rotational 

partition  function  (Q„)  is  very  temperature  dependent  and  is  given 

K 

by  [McClatchey,  et  al.  ] 


(6.2-2) 


where  the  exponent  Cj  assumes  values  dependent  upon  gas  type.  The 
vibrational  partition  function  (Q^)  is  also  temperature  dependent. 


but  not  in  such  a clear  manner.  Table  6.2-1  gives  the  value  of  C 


J 
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and  the  vibrational  functions  that  are  used  in  this  work.  For 


computer  implementation  it  is  simple  enough  to  use  the  different 
values  for  Cj  from  a look-up  array,  however  a speedy  and  accurate 
method  is  needed  for  the  vibrational  partition  function. 

Therefore  regression  analysis  of  the  vibrational  points  was 
performed  for  linear,  quadratic,  and  cubic  fits  to  determine  the 
best  relation.  For  each  gas  the  following  functions  of  temperature 
(T)  were  tried. 


Linear... 0^  = A^T  + A2 
Quadratic. . .Q^  = + B2T  + B^ 

Cubic... + 0.^'^  + C^T  + 

Linear  fits  were  abandoned  after  the  partition  function  was  graphed 
and  several  lines  were  tried.  For  the  quadratic  and  cubic  relations 
a computer  routine  was  written  to  determine  the  coefficients  B^  and 

^1- 

In  matrix  form  the  vibrational  partition  functions  can  be  written 
Q = AT,  (6.2-3) 

V 

where  is  a (7x1)  partition  array,  A is  a (7  x 3)  unknown  array, 
and  T is  a (3x1)  temperature  array.  For  cubic  fit  A is  (7  x 4)  and 
T is  (4x1).  The  unknown  array  A is  given  by  [Burr^] 

A = [T*^T]”^  (6.2-4) 
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where  the  superscript  (t)  indicates  matrix  transpose.  The  results 
of  the  regression  analysis  are  given  in  Table  6.2-2. 

In  view  of  the  small  errors  associated  with  the  quadratic 
coefficients,  they  are  used  in  the  transmittance  program  for  the 
vibrational  partition  function. 


6.3  Euler's  Solution  of  the  Transmittance  Differential  Equation 
The  transmittance  differential  equation  may  be  written 

^=crT,  (6.3-1) 

where 


a 


EK.M^ . 

J J J 


If  the  atmosphere  is  stratified  and  every  layer  has  a different 
o,  the  transmittance  relation  becomes 

t 


^1- 

dP  ‘^il+l  ’^J.+l' 


(6.3-2) 


where  dP  is  the  pressure  difference  between  the  i,  and  i+1  layers. 

This  is  a first-order  ordinary  differential  equation.  A solution, 

T (P) , is  desired  that  satisfies  Eq.  6.3-2  and  a single  specified 
initial  condition.  In  general  it  is  impossible  to  determine  t (P)  in 
exact  analytical  form.  Instead,  the  Interval  in  the  Independent 
variable,  P,  over  which  a solution  is  desired,  [0,P^],  is  divided  into 
sub-intervals  or  steps.  These  steps  correspond  to  the  pressure  differ- 
ences dP. 

An  application  of  the  Runge-Kutta  algorithms  to  Eq.  6.3-2  results 
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Molecule 


-9.88-10 

-2.16-10' 

-2.54-10" 

-2.81-10" 

0.0 

-6.00-10" 

-9.89-10" 


Quadratic  Coefficients 


4.49- 10 
6.77-10" 
5.16-10" 

4.28- 10" 

0.0 

2.28- 10" 

4.49- 10" 


Error  of 
Quadratic 
Formula 


9.6-10”^% 


2.1-10"^? 


1.3-10"^Z 


2.5-10"^?: 


0.0  Z 


i.3-io“^s: 


9.6-IO'^J: 


Molecule 


-1.68-10 


1.08-10 


-1.67-10 


8.85-10 


-6.38-10 


-6.91-10 


-1.68-10 


1.16-10 


-2.98-10 


1.16-10 


-9.46-10 


4.77-10 


-8.19-10 


1.16-10 


-2.63-10 


8.76-10 


-2.63-10 


2.05-10 


-1.16-10 


1.01-10 


-2.63-10 


Error  of 
Cubic 
Formula 

1.02 

8.9-io“^s: 

0.97 

2.2-lO'h 

1.02 

3.9  X 

0.99 

2.9-10"^X 

1.00 

2.8-10"®X 

0.99 

1.3-10"^X 

1.01 

8.9-10”^X 

Table  6.2-2.  Regression  coefficients  for  quadratic  and  cubic 
fit  of  Vibrational  Partition  Function. 
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In  the  improved  Euler's  method  or  Heun's  method  f Carnahan,  et  al.  ] 
of  solution  of  the  differential  equation.  The  desired  solution  of 
Eq.  6.3-2  is  thus  given  as 


dP 


(6.3-3) 


where 


^£+1  " "£  + • 


(6.3-4) 


subject  to  the  initial  condition  (top  of  the  atmosphere) 


t(o)  = 1 . 


(6.3-5) 


Euler's  algorithm  may  be  viewed  as  a predicting  equation  for 


-1 


(the  first  approximation  to  whereas  Eq.  6.3-3  may  be  con- 


sidered a correcting  equation  to  produce  an  improved  estimate  of  Tj^+1, 


Equation  6.3-3  is  used  iteratively  to  produce  a sequence  of  corrected 


values.  A maximum  of  ten  iterations  per  atmospheric  layer  is 


allowed  in  the  program. 

This  procedure  results  in  values  of  transmittance  for  every 
layer  of  the  atmosphere,  beginning  at  P=0,  t»1  and  continuing  to  the 
bottom  of  the  atmosphere  (ground  level) . The  variable  a must  of 
course  be  computed  for  each  layer  of  the  atmosphere  utilizing  the 
proper  equations.  However,  the  calculation  of  a is  sovereign  of  the 
solution  T,  and  is  done  independently. 


6.4  Transmittance  Iteration  Technique 

This  section  explains  a technique  that  was  developed  to  enable 
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the  transmittance  program  to  calculate  either  a monochromatic 
transmission  or  an  averaged  channel  transmission.  In  addition, 
the  option  to  extend  the  channel  limits  was  Incorporated  Into  the 


program. 

In  calculating  the  transmission  at  a single  frequency  the 
effects  of  all  lines  within  a specified  range  of  that  frequency  are 
added  together.  This  range  Is  denoted  by  the  variable  SEARCH.  When 
calculating  the  transmittance  of  a channel,  the  distance  BOUND  is 
used  to  extend  the  limits  of  the  channel  and  thus  account  for  lines 
that  contribute  to  absorption  within  the  channel,  but  whose  center 
frequency  Is  out  of  the  channel  (see  Fig.  3.4-3).  To  calculate 
the  transmittance  for  a band  the  monochromatic  transmit tances  are 
calculated  at  several  frequencies  and  then  averaged  over  the  band. 

The  following  variables  are  defined  below  and  used  in  the 
transmittance  program,  iteration  loop. 

WAVLO...  Lower  wavenumber  of  channel. 

WAVUP...  Upper  wavenumber  of  channel. 

SEARCH...  The  range  about  a single  frequency  within  which  the 
effects  of  spectral  lines  are  added  together. 

CNTR...  The  single  frequency  at  which  transmittance  is 
calculated. 

BOUND...  Amount  by  which  the  limits  of  a channel  are  enlarged 
to  consider  lines  that  absorb  within  the  channel,  but 
whose  center  frequency,  v^,  is  outside  of  the  channel. 

SPACE...  The  distance  between  the  single  frequencies  at  which 
transmittances  are  calculated. 
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INTVLS...  The  number  of  intervals  the  channel  is  divided  into, 
minus  one. 

INCR...  The  number  of  intervals  the  channel  is  divided  into. 

Also,  the  number  of  transmittance  calculations  made 
within  a channel. 

The  increment  loop  proceeds  as  follows  during  the  transmittance 
program  execution: 

(1)  The  variables  BOUND  and  INTVLS  are  read  as  input  data. 

(2)  The  lower  limit,  WAVLO,  and  the  upper  limit,  WAVUP,  of  the 
channel  are  read  as  input  data. 

(3)  The  limits  of  the  channel  are  adjusted  by  the  amount  BOUND. 

(4)  The  absorption  line  parameters  (v^,S,E,a)  of  lines  between  the 
adjusted  WAVLO  and  WAVUP  limits  are  transferred  from  tape  to 
computer  memory. 

(5)  The  variable  SEARCH  is  set  equal  to  BOUND.  If  BOUND  was  zero, 
SEARCH  is  set  equal  to  half  the  channel  width. 

(6)  The  channel  is  divided  into  the  desired  number  of  intervals 

(INTVLS+1)  and  the  variable  SPACE  is  set  equal  to  the  distance 

WAVUP -WAVLO 

between  transmission  calculation  points,  SPACE  = — i'ntVLS ’ 

(7)  The  first  frequency  of  transmission  calculation,  CNTR,  is  set 
to  WAVLO.  If  only  one  transmittance  calculation  is  to  be  made 
in  a channel,  CNTR  is  set  to  the  channel  center  frequency. 

(8)  The  upper  limit,  UP,  and  lower  limit,  DOWN,  are  set  to  CNTR  + 
SEARCH  respectively.  This  causes  the  transmittance  calculations 


to  Include  lines  within  the  distance  SEARCH  of  CNTR. 
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(9)  Transmission  computations  are  made  for  the  frequency  CNTR. 

(10)  If  another  transmittance  calculation  is  to  be  performed  in  this 
channel,  the  value  of  CNTR  is  changed  to  CNTR  + SPACE  (go  to 
step  9)  if  not,  another  channel  is  done  (go  to  step  2). 


(11)  After  all  channels  are  completed  the  program  terminates. 
Figures  6.4-1,  6.4-2,  and  6.4-3  graphically  display  this  procedure 
for  three  increment  options. 


INCR  = 1 


INTVLS  = 0 
BOUND  = 0 

C3JTR  = (WAVUP  + WAVLO)/2 
SEARCH  = (WAVUP-WAVLO)/2 


CNTR 


SEARCH- 


frequency 


WAVLO 

(DOWN) 


WAVUP 

(UP) 


channel 


Figure  6.4-1.  A single  transmission  calculation  per 
channel.  Channel  limits  not  adjusted. 
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INCR  = 1 


INTVLS  = 0 
BOUND  = 5 

CNTR  “ (WAVUP  + WAVLO)/2 
SEARCH  = BOUND 


CNTR 


DOWN 


UP 


WAVLO’ 

1 

SEARCH 

f— 

WAVUP' 

^ BOUND 

>1 

1 

* > 

1 

N 

f 

1 

WAVLO 

4 — 


channel 


WAVUP 


frequency 


Figure  6.4-2.  A single  transmission  calculation  per  channel. 
Channel  limits  adjusted  by  BOUND. 
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INTVLS  =•  2 


INCR  - 3 
BOUND  = 5 

CNTR  ■ WAVLO,  WAVLO  + SPACE,  WAVUP 
SEARCH  *=  BOUND 

SPACE  » (WAVUP-WAVLO) /INTVLS 


CNTR, 


DOWN, 


CNTR, 


DOWN 


1 

SEARCH 


4, 


T 


UP, 


SEARCH 


CNTR, 


DOWN 


3 

SEARCH 


WAVLO' 


WAVUP' 


i SPACE > < SPACE 


WAVLO 
► BOUNDS 


WAVUP 
>B0UND4 


Channel- 


Figure  6.4-3.  Three  transmittance  calculations  per  channel. 
Channel  limits  adjusted  by  BOUND. 


-frequency 


6.5  Transmittance  Program  Flovchart  and  Utilization 


The  transmittance  program  written  as  the  final  product  of  this 
work  la  Hated  in  Appendix  E.  The  equations  and  relations  necessary 
to  calculate  the  llne-by-llne  transmittance  were  implemented  in  the 
most  efficient  manner  possible.  There  are  five  main  portions  of  the 
program.  First,  all  Initializations  and  one-time  computations  are 
made.  Second  the  calculations  and  data  transfer  peculiar  to  a chan- 
nel are  completed.  The  third  step  Is  the  Increment  loop.  It  sets 
how  many  and  where  (frequencies)  transmittance  calculations  are  to 
be  performed.  Fourth,  the  equations  that  are  a function  of  atmos- 
pheric level  (altitude)  are  processed.  The  fifth  and  final  step  Is 
the  line  loop.  It  considers  only  the  appropriate  spectral  lines  in 
the  transmission  calculations.  The  output  is  a set  of  transmittance 
values  for  each  layer  of  the  atmosphere.  A general  flowchart  is 
given  in  Fig.  6.5-1. 

A major  portion  of  the  program  transfers  data  or  determines 
which  spectral  lines  to  use  in  calculations.  The  absorption  coeffi- 
cient computation  occurs  in  the  line  loop  and  utilizes  the  equations 
indicated  in  Table  6.5-1. 

The  program  was  written  in  FORTRAN  IV  for  an  IBM  360/65  computer 
[20,  21].  It  utilizes  two  subroutines;  DATSET  to  transfer  data  from 
tape  to  core  and  ZVOIGT  to  compute  the  volgt  line  shape.  Both  sub- 
routines are  explained  In  other  sections  and  have  listings  in  the 
appendices.  The  job  control  language  (JCL)  necessary  to  read  the 
data  tape  is  listed  In  Appendix  D with  the  DATSET  routine. 

The  card  Input  required  to  execute  the  program  Is: 
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NCHNL... 

Number  of  channels 

LEVS . . . 

Number  of  atmospheric  layers 

BOUND. . . 

Extentlon  amount  of  channel  limits 

INTVLS . . . 

Number  of  channel  intervals. 

minus 

one 

PRESS (LEVS)... 

Pressure  for  each  layer 

TEMP (LEVS) . . . 

Temperature  for  each  layer 

H20MIX(LEVS) . . . 

Water  vapor  concentration  for 

each 

layer 

03MIX(LEVS) . . . 

Ozone  concentration  for  each 

layer 

WAVLO... 

Lower  channel  limit  (for  each 

channel) 

WAVUP. . . 

Upper  channel  limit  (for  each  channel) 

Figure  6.5-2  shows  the  input  format  and  order  required. 

The  transmittance  program  was  executed  for  the  following  input 
data.  NCHNL=7,  LEVS=33,  BOUND«0,  INTVLS=0 

The  channel  limits  for  each  channel  were  taken  from  Table  3.4-1. 

The  pressure,  temperature,  water  vapor  concentration  and  ozone  con- 

32 

centratlon  were  taken  from  McClatchey,  et  al.  and  are  similar  to 
Table  5.3-3.  The  results  of  this  monochromatic  layer  transmittance 
for  each  channel  are  given  in  Appendix  G.  Other  computations  were 
made  utilizing  several  frequencies  per  channel  (INTVLS  ^ 0)  to  test 
the  Iteration  portion  of  the  program.  The  transmittance  of  a single 
channel  using  INTVLS  = 10  is  given  in  Appendix  H for  comparison  with 
the  monochromatic  calculation. 
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Card  Name 


Format 


Number  Required 


NCHNL 

1 1 1 

LEVS 

1 1 1 

BOUND  INTVLS 

Level-Channel  Data  1 | | 

LM 

M 1 

IX  13  IX 

13  IX  F10.3  IX  14 

PRESS(I)  TEM(I)  I 

II  I II 

i20MIX(I)  03MIX(I) 

Atmosphere  Data 

1 1 LEVS 

6X  E9 . 3 F5 , 

,1  IX 

E9.3  IX  E9.3 

WAVLO 

WAVUP 

Channel  limits  Data 

1 NCHNL 

IX  F10.3  IX  F10.3 

I 

[ 


Figure  6.5'-2.  Input  data  format  and  number  required. 
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CHAPTER  VII 

DISCUSSION  AND  CONCLUSION 

7.1  Preface 

This  final  chapter  briefly  comments  on  this  research  and  draws 
some  conclusions.  It  was  the  object  of  this  work  to  develop  a fast 
monochromatic  transmittance  program.  That  has  been  done.  The 
chapter  ends  with  a note  on  what  future,  related  works  might 
investigate. 

7.2  Discussion 

This  work  was  approached  with  the  intent  of  developing  a line- 
by-line  transmission  model.  To  construct  that  model  the  physics  of 
the  atmosphere  and  some  molecular  theory  of  gases  was  introduced. 

It  is  hoped  that  the  development  of  the  transmittance  algorithm 
proceeded  in  a fluid  manner  and  with  enough  detail  to  enable  the 
reader  to  follow. 

During  the  development  of  the  equations  that  were  needed,  an 
algorithm  to  quickly  evaluate  the  Voigt  function  was  born.  It  is 
sincerely  hoped  that  this  Voigt  algorithm  proves  to  be  a valuable 
one  for  atmospheric  work. 

7.3  Conclusions 

The  transmittance  program  written  as  a result  of  this  work  was 
used  to  compute  the  monochromatic  transmission  for  channels  in  both 
the  15  and  4.3  micrometer  band.  Also,  averaged  transmittances  were 
computed  for  selected  channels.  In  all  cases  the  program  executed 
properly  and  returned  values  that  seemed  to  be  correct. 
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The  results  of  the  program  were  not  compared  with  any  other, 
as  they  were  not  available.  However,  the  extinction  coefficients 
computed  in  the  program  were  compared  to  those  computed  in  a 

40 

completely  different  manner;  in  a program  written  by  Mr.  T.  Potter. 

The  agreement  of  the  coefficients  computed  in  the  two  different 
methods  was  excellent.  From  these  coefficients  it  is  a simple 
matter  to  compute  transmittance,  so  it  is  assumed  that  if  the 
extinction  coefficients  were  correct  the  transmit tances  would  be 
correct. 

The  single  most  important  conclusion  is  that  the  Voigt 
algorithm  developed  in  this  work  is  a rapid  and  accurate  method  to 
evaluate  the  Voigt  profile.  The  results  of  tests  that  were  made 
with  the  algorithm  surpassed  all  expectations  (see  chapter  four). 

7.4  Ongoing  Research 

It  is  my  opinion  that  the  Voigt  algorithm  developed  here  can  be 
Improved.  First,  the  routine  need  not  be  double  precision.  I be- 
lieve that  it  is  possible  to  obtain  very  good  results  with  a single 
precision  routine.  Secondly,  the  series  used  for  the  algorithm  in 
Region  I could  be  expanded  to  thirty  terms  and  a Chebyshev  econo- 
mization used  to  reduce  the  number  of  terms.  This  would  increase 
speed  and  maintain  the  desired  accuracy. 

The  results  of  this  transmittance  program  need  to  be  compared 
with  some  other  results.  Although  the  values  appear  to  be  correct, 
they  need  confirmation.  Once  the  validity  of  the  method  is  established, 
additional  data  needs  to  be  obtained  so  the  transmittance  of  the 


133 


atmosphere  can  be  more  closely  tabulated.  For  this  the  atmosphere 
should  be  divided  into  100  layers.  Additionally,  the  transmittance 
computed  for  a channel  needs  to  be  convolved  with  a slit  function 
corresponding  to  the  HIRS  instrument. 

Finally  a more  efficient  routine  needs  to  be  developed  to 
yield  an  averaged  transmittance.  The  method  used  in  this  work  is 
very  time  consuming.  To  calculate  a single  monochromatic  trans- 
mittance takes  about  15  CPU  seconds.  If  a step  of  0.01  wavenumbers 
is  used  for  a channel  with  a width  of  20  wavenumbers,  it  means 
2000  single  evaluations  or  about  8 hours  of  computation  time.  A 
better  system  must  be  developed. 
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APPENDIX  A 


KIELKOPF'S  VOIGT  ALGORITHM 


The  FORTRAN  IV  subroutine  that  implements  Kielkopf's  Voigt 
algorithm  is  listed  on  the  following  page.  It  is  a double  preci- 
sion routine  that  has  three  parameter  inputs  and  a single  parameter 
output.  The  input  parameters  are 


BL  = a. 


BG  = ajj/(ln2)  ' 

A V /I  9^1/2 

and  V = (ln2) 

“d 


The  value  of  the  Voigt  Function  is  returned  as  the  variable 
VGT  when  the  main  program  contains  the  following  statement 


CALL  C0EFF(BL,  BG,  V,  VGT). 


FORTRAN  tv  G LEVEL  21 


COEF 


DATE  - 77108 


11/01/30 


0001 

0002 

0003 

OOOA 

0005 

0006 
0007 
OCCA 

0009 

0010 
0011 
0012 
0013 
OOU 

0015 

0016 
0017 
0013 


SUBROUTINE  CCFr(3t.,DG,V,VCT) 

IM3|.i:tT  ReAL'‘f<(  A-U,0-/) 

DATA  rwO/2 .00100000/ 

A=HL/3G 

AL-»2./(  I . f . C9O*0L0G(  TWO)  *0  SORT  < (1.-.  090*0  LOG  I TWO) ) 

I ( (A  .'ODLCGt  TWO)  )/ « A*A)  ) ) ) 

G2  = t 1. /OLDOt  TWO)  ) /I l.-AL) 

8=HL/AL 

ETA=AL/ ( AL*G2) 

X=V/B 

GXsOEXP  t-(OLOG(TWO)*X*X ) ) 

ALXsl./d.tX^X) 

EX=  ( .80  29-.  A20  7*X*X)  /(  l.  + .203  0*X*X+.  0 73  3 5>'X*X*X*X) 

UX= (I .-ET A)*GX»PTA*ALX+ET A*( I .-ETA ) *E X*t GX-ALX) 

T*l./t  1.4-.A70A7*A) 

X l-A*( .6l6a6*T-. 1699A*T*T*  1.  32 5 5A* T* T* T)  /( 3. I A 1 592654*SL ) 
VGT«XIi'UX*(  1.772A533  51/  ( BL/OG)  ) 

RETURN 

END 
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APPENDIX  B 


DRAYSON's  VOIGT  ALGORITHM 

The  FORTRAN  IV  function  subprogram  that  uses  the  four  regions 
proposed  by  Drayson  Is  a single  precision  routine.  The  variable 
VOIGT  Is  set  equal  to  the  Voigt  function  value  In  the  subprogram 
when  the  main  routine  executes  a call  statement,  l.e., 

SHAPE  - VOIGT  (X,Y). 


The  parameters  X and  Y are: 


Y - — (ln2)^^^ 


“d 

v-v 


X 


The  program  listing  begins  on  the  next  page. 
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FO^T!»A'<  IV  G tEVPt  21 


VOIGT 


OATF 


77108 


ll/Ol/JO 


OOOl  FUNCnON  vniGTIX.Y) 

03  )2  RF  AL  ‘^(22  )70..  .70T<602F-7/  , P I ( I 5 ) i XN(  I 5)  / 1 0.  , 9.  , 2*  8.  , 7.  , 6.  1 5 . , A.  , 

1 7<-3.7tYNIl5)/3*.6,.5,2*.A,A'>.3,l.,.9,.a»2*.7/, 

2 Dn(25),Ol(25),r)2l?5)in3l25),OA(2a),HM(25), 

3 H/  .701  / ,XX  (3  I / . 32A6A76  . I .65J69f  .7071  U6.87, 

3 HH(  3)/.  2562121, .25n-(26fiE-l,.28209A8/, 

A N 8Y  ?(  19)7  ^#.5, 9 »,8.5,  8.,  7 .5,  7«,  6.5,6.  ,5. 5, 5., A.  5,  A.  ,3.5,3.  ,2.5, 

5 2.  , 1 . 5, 1 . , .•>7, 

6 0(21)  7.  7093  6 02--7,-.2  5l8A3AE-6,.  35668 7A6-6  ,-.2  78  763  8E-5, 

7 . 3660 7 A E-5,-. 2 555  55  16-4, . 722 8 775F-4, - . 1 93 343 16-3 , . A8995206-3 , 

7 -.  1 1 73 2 67F -2  ,.2  646 762  6-2  ,-.56  23190  6-2, .11  19601  6-1 , -.20349766-1 , 

P .3621  5736- 1,-. 5 85 lAl  126- 1, . 87 70816“- I,-.  121664,.  1 558 A 18 A , . 2 7 


C003 

LOGICAL  T9'J7.FALS6.7 

OOOA 

IFt  TRU)  GO  TO  104 

00)5 

TRIJ=.TRUE. 

0006 

DO  101  1=1,15 

0007 

101 

RI  ( I )=-I 72. 

00  18 

DO  103  1=1,25 

0309 

H'4(  I ) I-.5) 

0010 

Cn=4.’!=HN(  I ) 7HN(  I )725.-2. 

001 1 

on  102  J=2  ,21 

00  12 

102 

ei  o*-i)=co«8(  j)~3(  j-i)  + C(j) 

0013 

D0(  I ) = HN I I ) *( BI 22)-e( 2 1) )75. 

OOIA 

Did)  =l  .-2.*HN(  I )’:<r)C(I  ) 

00  15 

D2(  I ) = <HN(  I 1 *D  K n + Oni  1 ) ) 7Rt  (2) 

0316 

D3 I I )= ( HN(  I ) *D2  t 1 ) +01 ( n )7Rl (3) 

00  17 

103 

04( I ) =( HN(  I l*03t  I ) + 02 ( I ) ) 7RI 14) 

00  18 

104 

IF(X-5.)  105,112,112 

0019 

105 

IF(Y-1.)  110,110,106 

00  20 

106 

1F(  X.r.T.  1.  85*(  3.  6-Y)  ) GO  TO  112 

0021 

IF  lY  .LT .1 .45)  GO  TO  107 

C022 

1 =Y  + Y 

0023 

GO  TO  138 

0024 

107 

I = 1 1 . =>Y 

0025 

108 

J = X+XH.85 

0026 

M4X=XVIJ l^YNI I )+.46 

00  27 

MIN=.MIN0(  16  ,21-2»MAX  ) 

0023 

UU=Y 

0029 

VV  = X 

0030 

DO  109  J=MIN,19 

00  31 

17^=NBY2IJ)  7(UU*UU+VV»VV) 

0032 

,UU= Y+U*UU 

3033 

109 

'vv=X-U=i'VV 

00  34 

V3IGT='JU7(  UU*UU+VV*VV)  71.7  72454 

0035 

RETURN 

0034 

110 

Y2  =Y*Y 

00  37 

IFI X+Y.GE.5. ) GO  TO  113 

0033 

N = X7H 

C009 

0X=X-HN{ N+l ) 

00  40 

U=(.(  1 D4(  N+  1)  *0  X+0  3(N+  1)  ) *'DX+02  (N+l)  )*0X+D1  ( N + l ) ) *DX + 00  ( N+ 1 ) 

1)41 

V=1\-2.’>X9U 

004? 

VV=EXP(  Y2-X«X)«-CCSI2.*X*Y  )71 .12  8379-Y*V 

00  43 

'JU=-Y 

0044 

•1AX=5.+  (12  .5-X  )*.fl*Y 

0045 

DO  III  1=2, MAX, 2 

0046 

U=(X«V+U)7RI{ I ) 

004  7 

V=(X*U+V)7Rl( I+l ) 

0043 

UU=-UU* Y2 

FC0TR4N  IV  G LEVEL  21 


VOIGT 


DATE 


77108 


n/Ol/30 


oo'fO 
00  50 
OCSl 
00*:  ’ 
10S3 
00  54 
00  5 5 

00  5 > 
00  57 
00  58 
0050 
0060 
)OSl 

0052 

0053 

0054 


111  VV-VV*V*UIJ 

vn  ir,T=  I . I28370«vv 
RE  TURN 

112  Y2^Y<'r 

IT  (Y  .LT  .ll.-.6875*»X)  GO  TO  113 
U=X-XX( 3) 

V=X*XXl 3) 

VO  IOT=Y  * (HH(3  ) / ( Y2  »U<'U  ) ♦MH  ( 3 >/  ( Y2  •■V'OV  ) ) 


RETURN 

113  U=X-XX(  I I 

v=x+xx( ; ) 

UU=X~XX( 2) 

VV=X^XX( 2) 

VniGT=Y*(HH(l»/(Y2«-U*U)*HHIl  )/ ( Y2 +V*V  ) i-HH  1 2 )/ C Y2 +UU  *00  J fr 


I HM(  2) /(  Y2»VV*VV)  I 
RETURN 
ENO 
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PROPOSED  VOIGT  ALGORITHM 

The  Implementation  of  the  proposed  Voigt  algorithm  utilizing 
three  regions  of  the  XY  plane  Is  a FORTRAN  IV  function  subprogram. 
The  routine  sets  the  variable  ZVOIGT  equal  to  the  value  of  the  Voigt 
profile,  for  the  parameters 


Y - ^ (102)1/2 
“d 


and  X - (ln2)^^^ 


“d 


The  call  statement  required  In  the  main  program  Is 


SHAPE  = ZVOIGT  (X,Y). 
The  program  listing  Is  on  the  next  page. 
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(V  <S  LEVEL  21 


Z vnicT 


DATE  - 77108 


11/01/30 


E 

► 

E 


t 


( 

E 


0001 

0002 

0003 


OOOA 

0005 

0006 
0037 
0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 
0017 
0013 
0019 
00  20 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 


FUNCTION  ZVCIGTIXiY) 

INPLICl T 3EAL*8( A-H,0-Z) 

KEAL*9  AN! 30)/ 1 . 0 OO ) 7000. 3 3333333» . I OOOOOOOt -.238095238E-1 , 

1 4.629  62963F-3,-7.5757575  7E-4,  I .068  376069E*>4,-1 . 322751323E«5t 

2 1.  45a9l69lE-6,-l.4503  8 522.'»';-7,l.  ;i22  53296E-8.-)  .089222104E-9, 

3 9.350702705E-Ut  -5.S477'’40:4'--nf  3.951i429516E-l3t 

4 -2.46682701^-14, 1.443326465E-1 5,-8.0327350128-1 7, 

5 4.22l407289E-in,-2. 10785 5 192E- 19, 

4 1.002516493E-20.-4.551046759E-22, 1.977064754E-23, 

♦-8. 230l49299E-25,3.239260349F-26,-l.264l07899E-27,4.67848352E-29, 
♦ -1.669761793E-30, 5.7541916448-32,-1.9169428628-33/, 

*■  Al/. 46131350/, 

6 A2/. 19016350/,A3/.C9999216/,A4/1. 78449270/, A5/.002883894/, 

7 A6/5. 52534370/, 81/ .5 1242424/ ,92/ . 27525510/, 83/. 05176536/ , 

8 84/2. 72474500/, PI  SO/1.128379167/ 

S=»X*X-Y*Y 

T-2.*X*Y 

IF(Y.GE.5..0P.X.G8.S.)  GO  TO  111 

IFIY.GE.l.S  .OR.  X.GE.3.0)  GO  TO  112 

XSER*Y 

YSER—X 

XN«Y 

YN— X 

X2--S 

Y2’— T 

N-6. 8424X43. 0 
IF(N.GT.29)  N-29 
IF(X.EO.O.O)  N-15 
DO  10  I 10-1,N 
XNEW«XN*X2-YN*Y2 
YNEW»Y24XN4YN*X2 
XSER-XSER4XNEW4AN(I 1041) 

YSFR*YSER4YNEV<*AN(  lio+l  » 

XN>XNEW 
YN»YNeW 
13  CONTINUE 

ZVOIGr»OEXPI-S)4lOCOS(-T)*Il.-PISO*XSER»4PISQ40SINI-T)*Y$E«l 
RETURN 
112  R-T4T 
T5T*X 
•F-S-A6 
• G-S-A4 
H-S-A2 

ZVOIGT«Al*l  IT-H4Y)/(H4H4R))4A3*I  IT-G*Y)  /IG4G4R))4> 

I A5*HT-F*Y)/ (F*F4R)) 

RETURN 
111  R«T*T 
T«T*X 
F-S-82 
G*S484 

ZVOIGT-81*! IT-F*Y)/IF4F4R) ) 483* I IT-GFY 1/ |G«G4R)| 

RETURN 

END 
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APPENDIX  D 


DATSET  SUBROUTINE 


The  data  processing  routine  used  by  the  transmittance  program 
to  transfer  data  from  tape  to  computer  memory  or  disk  is  listed  on 
the  next  page.  The  parameters  of  the  FORTRAN  IV  subroutine  are 
explained  in  section  5.6.  The  following  JCL  must  be  included  for 
program  execution.  It  defines  the  seven  files  of  the  tape  so  that 
the  routine  can  more  quickly  select  that  portion  of  the  tape  from 
which  to  read  data. 


//GO.FTllFOOl  DD  UNIT=2400,DSN=F1L1 , 

//  LABEL=(02,BLP,IN) ,DISP=0LD,V0L=SER=T00388 , 

//  DCB=(RECFM=F,BLKSIZE=3900) 

Seven  of  these  cards  must  be  used,  the  remaining  six  are  exactly  the 
same  except  as  shown. 

...  GO . FT12F001 . . . DSN=F1L12 . . . LABEL= (05 , BLP , IN) . . . 

. . . GO . FTl 3F001 . . . DSN=F1L1 3 . . . LABEL= (08 , BLP , IN) . . . 

...  GO . FT14F001 . . . DSN=F1L14 . . . LABEL= (11 , BLP , IN) . . . 

...  GO . FT15F001 . . . DSN=F1L15 . . . LABEL= (14 , BLP , IN) . . . 

...  GO . FTl 6F001 . . . DSN=F1L16 . . . LABEL- (1 7 , BLP , IN) . . . 

.. .G0.FT17F001. . .DSN-F1L17.. .LABEL=(20,BLP,IN). . . 

The  routine  can  transfer  3000  lines  of  data  to  core.  If  disk  Is  to 
be  used  the  necessary  JCL  and  DEFINE  FILE  statements  must  be  Included, 
and  in  the  DATSET  routine  the  comment  card  used  and  the 

SPECS (...)= 

and  KIND  (...)-  Cards 

removed. 
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F0PT3AN  IV  G LEVEL  21 


OATSET 


DATE 


77108 


11/01/30 


00  01 

5U300UTINE  nATS'=T(WAVLrj,WAVUP,  ICNT.MAX.SPECS.K  INO,  LTOTAL  ) 

0002 

01  MENS  ION  IDATFl  40)  , ! SOI  43)  ,'3(iL(  43)  ,G8G(  23)  , ICNT(7) 

0003 

01  Mr  NS  I ON  WAV  (40  ) ,STRC.H(40  ) , WOT  H ( 40  ) , ENGY(  40  ) . ROT  1 0«  40,5)  , VIB(  401 

OOOA 

0IM5NSI0N  SPEC S( 3030 ,4) ,KINO(3000) 

)0  )5 

LTOTAL=0 

OOOS 

DO  40  140=1 ,7 

0007 

40 

ICNT(  I 40)  = 0 

0008 

NF  IL=11 

0000 

(F{ WAVLO.GT. 500.0)  NF[L=12 

00  10 

IF( WAVLO.OT.  lOOO.O)  NFIL=13 

0011 

IF  (WAVLO.vST  .2003  .0  ) NFIL=14 

0012 

IF ( WA VLO.G T. 5000. 0)  NFIL=15 

00  13 

IF(WAVLO.GT .7503.3)  NF[L=16 

001  A 

IF( WAVLO.GT. 10000.0)  NFIL=17 

0315 

07 

REAO( NF IL ,01,EN3=99)  NUM.FREQ 

OOIG 

31 

format  ( I 1 3, FIO .3 ) 

0017 

IF ( FREQ. LT. WAVLC)  GO  TO  07 

0013 

backspace  NFIL 

0310 

BACKSPACE  NFIL 

0070 

23 

PFAO(  MF  I L ,1 1 ,FNi3  = 9 3 ) NUM  ,(WAVtK),STRGH(K),WOTH{K), 

1 FNGY(K),<R0TI0(K,J),J=1,5),VIB(K),I0ATE(K) ,ISO(K), 

1 VOL  (K  ) , K = 1 ,40  ),  (GtlGlM  ) ,M=1 ,23  ) 

0021 

11 

FORMAT!  I 13,40(F10.3 ,E10.3 ,F5.3,F10,3,5A6, A5,  13,  14,  I3),23A30) 

0022 

on  10  113=1, NUM 

C023 

KP0L  = M0L(  I 10) 

3024 

MOLF=  lO^KMOL 

0025 

IF(WAV(  113)  .LT  .WAVLO)  GO  TO  10 

0026 

IF ( WAV(  I 10)  .GT.WAVUP)  GO  TO  13 

0027 

ICNTIKMQL )=ICNT( KMnL)+l 

C 

WRITE (MOL  F*  ICNT (KMOL ) ) WAV! I 10 ) ,STRGH!  I 10) , WOTH! 1 10) ,ENGY(  110) 

0023 

L T0T4L  = L TOTAL* 1 

0020 

IF (LTOTAL  .GT .MAX ) GO  TO  97 

0030 

SOECS! LTCT AL,1 ) =W AV ( 1 10) 

0031 

SPEC  S(LT0*AL,2)=STRGH!I  10) 

0032 

SPrCS  (LTOTAL,  3 l = W3rH(  110) 

3033 

SPECS! L TOTAL ,4 ) =fnGY ! 110) 

00  34 

KIND!LT0TAL )=M0L! 110) 

0035 

1 0 

CONT  INUE 

00  3 6 

G.3  TO  2 3 

0037 

13 

continue 

00  33 

Rewind  nfil 

0C39 

• 

, GO  TO  999 

004C 

99 

WR  1tF(6, 21 ) NF IL 

0041 

21 

FORMAT ! 1 X ,35HLOWER  WAVE  NUMBER  NOT  FOUND  ON  TAPE,/, 

1 1X,15HFIL6  number  I S: ,I10,/,1X,18HPROGRAM  TERMINATED,//) 

0342 

STOP 

00^3 

98 

WRITE!6,31)  NF I L,WAV (40) ,W AVUP 

0044 

31 

FORMAT! //, 2X, 31MFN3  OF  FILE  ENCOUNTERED  ON  F I LE , 1 3 ,/ ,2X , 

1 34HBFF0qF  UPPER  WAVE  NUMBER  WAS  FOUND, /,2X, 

2 ^OHLAST  wave  NUMTfa  I S:  ,F10.3,3X, 

3 25HUPPER  LIMIT  S^FCIFIFO  WAS , 5X, F 1 0. 3 , //2 X , 

4 18HPR0GRAM  CCNT INUING,/// ) 

0045 

REW  IND  NF  IL 

3346 

NF IL=NF  tL*l 

0047 

IF!NFIL.GE.17)  NFIL=17 

0048 

GO  TO  23 

C049 

97 

wai TE (6  ,4  1)  1 I CNT !M) ,M  = 1 , 7) ,WAV!40) 

0050 

41 

FORMAT! //,  IX, 31HE XPECTEO  PECORO  NUMBER  EXCE EOEO,/ , IX , 

1A3 


pr.pT^AN  IV  r,  LEVFL  21  OATSET  DATE  = 7710d  11/01/30 


1 3-^HNUMBE‘l  OF  BECQPDS  OF  EACH  TYPE  I S : , / ,2 X ,7  ( 16  .2X  I , / . 

2 EOHLAST  WAVE  NUM'Jcq  | S : f 2X  , F 1 0. 3t  / . 2 X,  1 8HPR06RAM  CONT I NU I NG  i ///) 
ICNTI KMOL) »ICNT IKMOL)-l 

REWIND  NF I L 
999  RETURN 
END  . 


0051 

0032 

0353 

COS'* 


/ 


The  transmittance  program  written  aa  a result  of  this  work 
follows.  The  Important  variables  of  the  routine  are  defined  In 
Appendix  F.  The  logic  and  a flow  chart  for  the  program  Is  contained 
In  Chapter  six. 

f 

i 

i 


i 

1 


FOPTrAN  IV  '■i  UFVEL  21 


MAIN 


DATE 


77108 


11/01/30 


0001 

IMPLICIT  REAL’i'n  ( A-n,r:-2) 

0002 

REAL  WA  VLO.ViA  VUP  ,F  RED, NIGH,  SIZE  ,EGY,  ROUND 

0003 

c 

RFAL>!-6  TE'«P(033)  , PRFSS  (033>,CJt7)/l.5,l.0,l.5,2*l.0,1.5,l.0/, 

1 XT  ^OLI  7) / 13.  , 66.  ,68.  ,66.  ,?  i.  , 1 6 . , 32  . / , N? OM IX ( 03 3 ) . 05  MI  X 1 0 33 ) 
1 SPECS  (300  0,6  ) 

00  OA 

I NrEr,ER»6  KI  NOt  lOOOl  ,MAX/3')00/ 

0003 

DIMENSION  0VI33, 7) ,A8SX( 33,71 ,ICNTI7)  , TRANS (33)  , X 1 I 33 » ,X 2 1 33 ) , 
1 X3(33) 

OOOS 

DATA  P 0/1  01 3.1  /I  2 2 6:10/  ,T  0/ 295. 00000000/ , HCKT  0/6 .860728  2 16E-0  3/ 

0007 

DATA  ALG2/ 1.83 255661 12/, SOHPI V/0. 5661895335/ 

0008 

READ(5,01)  \C  (\L, LPVS , aOUNO,  IMTVLS 

000  9 

01 

FORMAT!  IX,I3,IX,I3,1X,F10. 3,1X,I6) 

0010 

DO  11  110=1, LEVS 

00  ll 

10 

RE  AO ( 5, 1 1 ) PRE  SS( t 1 01  , TEMPI  I 10) ,F2CMIX { IlO I,03MIX(  I 10) 

00  12 

11 

FORMAT! 6X, F 9.3. F 5.1,lX,Fg.3,lX,E9. 3) 

C 

c 

CNE  TIME  CALCULATIONS 

0013 

L 

on  160  1160=1, LEVS 

00  1 A 

T2  = TEMPI  1160)=>TFMPni60) 

0015 

0V(  I 161, 1 ) = .99-316626+.669331030-06*TEM?(  I 1 60 ) -.  98896  03  OO-OT*  T2 

0016 

0V(  I 160  ,2)  =1 .0  7075l3*^.51691620n-03*TEMP(  1160  )- . 2568  366  lO-O  5#T  2 

0017 

Q VI 1 161, 3) =.5896  9668^ . 67  78/5030-03*  TEMPI  1 1 60 ) -. 2 1 57 3 75 1 0-05*T2 

1018 

OV  (I  15 1,61  = 1.1211 63  3 + . 6 282 30C0O-O3*TEMP(  1 160)-.  281 01 3 720-05* T2 

0019 

0V( I 160,5) =1  .01000000 

0020 

0V( I 160, 6)  = .98 56 3537+. 2 2 7665130-03* TEMP( I 160)-. 6 006993  30-06*12 

0021 

OV  (1 160,7 )= .99 5066 26 + .669 33 1030-06* TEMPI  I 150 )-, 9 389603 90-07*T2 

0022 

Xl(  ! 160)  = T0/TEMP(  I 160) 

0023 

X2(  115)  )= ( PRESS!  I 160)/P0) *OSORT( Xl(  I 1501  1 

002'* 

X3 ( 1 150) =HCKTO* ( 1 .-XI (1160) ) 

0023 

160 

CONTINUE 

U 

C 

c 

CHANNEL  LOOP 

00?3 

00  20  120=1, NCHNL 

0027 

B5An{5,2l)  WAVLO.WAVUP 

0023 

21 

FORMATI 2(  1X,F10.3)) 

0029 

WA VLO=WAVLO-BOUNO 

00  30 

WAVU°=WAVUP+aOUNO 

0031 

Seapch=  sound 

003  2 

IF(  30'JNO.EC.O.O)  SeARCH=  (W  AVUP-WAVL0)/2  .0 

0033 

CALL  OATSETt WAVLO.XAVUP ,ICNT ,MAX , SPE CS , KI NO , LEND) 

003A 

INCR=INTVLS+l 

0035 

IF( INTVLS.NE.O)  SPACF=(WAVUP-WAVLO-2.*BOUNO)/ INTVLS 

0036 

on  60  161=1, LEVS 

0037 

TRANS( I60)=0.0 

00  38 

c 

60 

CONTINUE 

C 

C 

r 

increment  LOOP 

0039 

Q 

IFIRST=1 

00%  0 

CNTR=WAVLO+QCUND 

0061 

IF( INTVLS.EC.O)  CNTR=(WAVL0+XAVUP)/2.0 

00A2 

DO  60  160=1,  INCR 

006  3 

UP=CNTR+SEARCH 

0066 

OOWN=CNTR-SEARCH 

0065 

DO  70  I70=IFIRST,LEN0 

PC’  T'^AN  IV  <i  LEVI’L  21 


MAIN 


l-A  r 


7/ioa 


1046 

CHfCK^SP  ;cs ( I 70, 1 ) 

00  4 7 

1 1' (c‘i‘ cs.on.DusN)  on  to  i3 

00<.8 

70 

r.d'i  r iNu*-: 

1149 

: '•  r < st 

00  50 

r.O  i(!  2 3 

1051 

13 

iEK5r=  17) 

0052 

23 

CON  I I NllE 

00  53 

DO  30  I 30= IP  JR  ST, LEND 

0054 

CHrrK=S“ECS  i MO,  1 ) - 

0055 

IF  (CHPC.S.GE.liP)  GO  TO  33 

0056 

80 

CONTINUE 

0)57 

ILAST=LEN0 

0053 

GO  10  43 

0059 

33 

II.  AST=  I 80 

0)6) 

43 

CONTINUE 

C 

LEVEL  LOOP 

0061 

C 

DO  30  I 30  = 1 ,LE VS 

0062 

00  50  159=1,7 

0063 

AOSKI  130,  I5i))=0.0 

0064 

#■ 

50 

CONTINUE 

Vtf 

C 

c 

LINE  LOOP 

005  5 

DO  90  190= IF IRST , ILA5T 

0065 

ITYPE=KI NOI 190  1 

00  57 

A0  = 3. 539-0  7*SOR7( TFMOII 301 /WTMOL ( 1 TYPE) ) 

0068 

X4-X  11  : 3)  ) ^^CJ  1 [TYPE)  <=QVI  130,  ITYPE) 

0069 

ALFAD  = AD»SPrCS(T90,l  ) 

0070 

ALF4L=X2(  130M'SPECSn90,3) 

0)71 

YFAC= ALG2/ALFAD 

0072 

U=ALrAL* YEAC 

0073 

RKO=  I')EXP(X3n30  )*5PECS(  190,  4)  » X=SPEC S C 1 90,  2 ) *'X4)*YFAC 

0074 

Y=0A3S(  CNTR-SPr  CSI 190,1 ) )»YFAC 

0075 

SHAPe=ZVOIGT( Y,U) 

0076 

ABSKinO,  ITYPE)  = ABSK(  130,  ITYPE  ) ♦RKO*SHAPE 

0077 

90 

CONTINUE 

0078 

c 

30 

continue 

c 

ABSORPTiON  CGEFFICEINT  ♦ MIXING  RATIO 

ccr9 

0080 


0081 


0082 

OCSJ 

0034 

0085 

0086 
0087 


C 

c 

c 


DO  130  1130=1, LEVS 

ABSKI  I I 30,  n = ( I 130,  U*2.7434  75Aqon;6*H2CMIXI  U30)«- 

1 ABSKC  Il31,2)='6.0546e2AA70l5tA8SK(  t 1 30,  3 ) * 1 . 279623310  lO-* 

2 03MIXI  n30)t-AaSX(  ^130,4)♦5.89  703  5853Dl^^■ 

3 ARSK  I I 130,5)  *1. 5894307  130  12t- AB  SK I 1 1 30, 6)  *3.  381698754013  «• 
A ABSK(  Il30,7)«4.443;3'?014016)*ISQRPlVi 

130  CONTINUE 

TR ANSM ITT ANCE 

LSTPP=LPVS-1 
TAU=1. 

00  150  I l50-->l,LSTOi> 

INUM=1  ^ 

HSTPR  = ( press  1 1 150*U  -PirSS  < i 150  ) )XlO00  . 

TPREO=TAU-"STtP’^  \ri  SK(  I If-O,  U =»  TAU 


./11/30 


F0=!TRAi'4  IV  G LEVEL  21  MAIN  DATE  = 77108  11/01/30 


0T38 

63 

TALI\'EW=  14U-(*iSTrP/2.1*(A(iSK{;i50tl  1 «TAU»  A3  S K I I 150  ♦ V t i • > PREOi 

OO:!*! 

IF  I DA -iS  ( T P'’r  )-r  A'JNEW  1 .LE.O  .0001  .QR  . INUM  .Gc  . 10)  GO  VO  53 

00«70 

IPREO-  rAUNEW 

0001 

CC*>2 

GO  m .jO 

0003 

53 

•IAU=  ■ AUNE  W 

O-OOA 

K-tl5)«^l 

0005 

Ti-  ANS  1 .<  1 =THA,\S  I K)  «-T  AU 

0C«6 

150 

CON  riNU^ 

0007 

CNTF^  CMRtSPACE 

0003 

50 

continue 

0090 

TRAN<;i  n = INCR 

01  00 

WAVl.  O=WAVI.'0+3CUNO 

OVOl 

WA  vi!t*=w\vup-aouNr) 

0102 

WRITE(6, 31)  1 ?0.WAVU0tW4VUP 

01  03 

31 

FORM  AT  1 1 HI ,5 IX . 1 AHCHANNEL  NUMBER, IX,  1 2, / , 23X, 

1 23H;t'^GtNNlNG  XA  V-NUMIER  OF  , 1 X , F 10 . 3 , 1 X , 3 H AND,  1 X, 

2 20HFr;oi'iG  WAVr‘i'J-1Br;«  OF  , lX.P  10.  3, /,  25X,  701  IH*  ) ,/,25X,  701 IH*)  ,/) 

01  JA 

WKilF (5  ,101  ) 

0105 

LEVEL=LFV5 

01  06 

DO  150  1150=1, LEVS 

0107 

TRANS! 1 150) =TRANS( 1 150) /INCR 

0106 

WR  IFE(6,91)  LEVcL.PRESSI 1 150) .TEMP  1 1 150)  .TRANS! 1150) 

0100 

LEVEL=LEVS-I150 

ono 

101 

FORMAT!  /,19X,5HLEVEL,UX,8H  PRES  SURE,  15X,llHT  EMPERATURE,  9X, 

1 22HAVERAGeD  TR AN SM 1 T TANC E , / , 1 6X, 85! IH= ) ) 

Clll 

91 

FORMAT!20X,I3,11X,010.5,17X,F5.1,15X,017.10) 

0112 

150 

CONTINUE 

0113 

20 

COM  INUE 

0115 

STOP 

0115 

END 

APPENDIX  F 

DEFINITION  OF  VARIABLES  USED  IN  THE  TRANSMITTANCE 

PROGRAM 

The  following  is  a list  of  the  important  variables  and  their 


definition. 

used  in  the  transmittance  program. 

ABSK(33,7) 

Array  of  absorption  coefficients 

ALFAL 

“l 

ALFAD 

“d 

A1G2 

(ln2)^/^ 

BOUND 

Amount  by  which  channel  limits  are 
changed . 

CJ(7) 

Array  of  rotational  partition  function 
exponents. 

CNTR 

Frequency  at  which  transmittance  is 
calculated. 

DOWN 

Lower  limit  of  the  range  of  lines 
included  in  transmittance  calculations. 

EGY 

Energy  of  a line  (E)  = SPECS(I,4). 

FREQ 

Line  wavenumber  (Vq)  ■ SPECS (1,1). 

HIGH 

Line  strength  (S)  ■ SPECS (1,2). 

HSTEP 

dP  in  bars. 

H20MIX(33) 

Array  of  water  vapor  concentrations. 

IFIRST 

Subscript  of  the  lower  limit  line  used 
in  transmittance  calculations. 

ILAST 

Subscript  of  the  upper  limit  line  used 
in  transmittance  calculations. 
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i 


! 


I 

I 

} 


1 

1 

1 


INTVLS 

Number  of  Intervals  In  a channel. 

KIND (3000) 

Array  of  molecular  types. 

LEVS 

Number  of  layers  of  the  atmosphere. 

NCHNL 

Number  of  channels. 

03MIX(33) 

Array  of  ozone  concentrations. 

PO 

Reference  pressure  = 1013. 17mb. 

PRESS (33) 

Array  of  layer  pressures. 

QV(33,7) 

Array  of  vibration  partition  functions. 

SEARCH 

Range  from  CNTR  within  which  effects  of 
lines  are  added  together. 

SHAPE 

Voigt  profile  value. 

SIZE 

Line  width  (a^)  = SPECS (1,3). 

SPACE 

Distance  between  points  where  transmittance 
is  computed. 

TAU 

Iterative  transmittance. 

TAUNEW 

Corrected  iterative  transmittance. 

TEMP (33) 

Array  of  layer  temperatures. 

TO 

Reference  temperature  = 296. 

TPRED 

Predicted  transmittance. 

TRANS (33) 

Transmit tances  of  each  layer. 

U 

“d 

UP 

upper  limit  of  the  range  of  lines  in- 
cluded in  transmittance  calculations. 

WAVLO 

Channel  lower  limit. 

WAVUP 

Channel  upper  limit. 

WTMOL(7) 

Molecular  weights  of  atmospheric  gases. 

A 
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APPENDIX  G 

■ S 

MONOCHROMATIC  TRANSMITTANCES 

The  results  of  monochromatic  transmittance  calculations  for 
the  seven  channels  of  the  HIRS  are  listed.  The  frequency  at  which 

■ t 

the  transmittance  was  computed  was  the  channel  center-frequency. 
Only  lines  between  the  channel  limits  were  used. 
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CHANNEL  NUMBER  I 

BEGINNING  WAVENUMBER  OF  665.200  ANO  ENDING  WAVENUMBER  OF  670.800 


■VEL 

prf:ssure 

TEMPERATURE 

AVERAGED  TRANSMITTANCE 

t:?  =s  3 

SSX3«S«3S'5SSS3SS-~XS 

aasa=a*a  = aa3s=--  = »i  's-  3»aia»aa3  = 

3XSS3XXSSSX3SX9SX3SXBXXSI 

33 

0.  3C00E-O3' 

2 1 0.  0 

).l OD )J0 IT JOO  31 

32 

0.5710E-01 

218.0 

0.9'<9999  .320  00 

31 

O.^-SIOE  00 

276.0 

0 .9  /9939  T9350  00 

30 

0.176  3E  01 

270.0 

0.99997533330  CO 

2<J 

0.3J30E  01 

258.0 

0.9999663063D  00 

28 

0.  65206  01 

265.0 

0.9993706895D  00 

27 

0.l320e  02 

236.0 

0.99970339910  00 

26 

0.2770E  02 

226.0 

0.99938633670  00 

25 

0.3220E  02 

223.  0 

0.99933258360  03 

2A 

0.3760E  02 

222.0 

0.99921257620  00 

23 

0.  4370E  02 

220.0 

0.99911961300  00 

22 

0.5133E  02 

219.0 

0.99901 778600  00 

21 

0.5950E  02 

218.0 

0.99893873500  00 

20 

0. 6950E  02 

217.0 

3.99379309600  03 

19 

3.8120E  02 

216.0 

0.99866066700  00 

18 

0.9500E  02 

216.0 

0.99351623960  00 

17 

O.lllOE  03 

216.  0 

0.99335706670  03 

16 

0.1300E  03 

216.0 

0.99817532630  00 

15 

0. 1530E  03 

216.0 

0.99796658620  00 

16 

3.17S0E  03 

216.  0 

0.99773605050  00 

13 

0.2090E  03 

222.0 

0.99767757100  00 

12 

0.  26  3 0E  03 

229.0 

0.99719837260  00 

1 1 

3.2310E  03 

235.0 

0.99689389680  00 

10 

0.3260E  03 

262.0 

0.996573o7980  00 

9 

0.3723E  03 

268.  0 

0.99622532110  03 

9 

0.6260E  03 

255.0 

0.99586771810  00 

7 

0.667CE  03 

261  .0 

0.99563699350  00 

6 

3.5560E  03 

267.0 

0.99530102660  00 

5 

0.&280E  03 

273.0 

0.99653662030  00 

6 

0.  7100E  03 

279.0 

0.99603218660  00 

3 

3.8020E  03 

285.0 

0.99368372730  00 

2 

0.9020E  03 

290.0 

0.99290212630  00 

1 

0. 1313E  06 

296.0 

0.99227078150  00 

CHANNEL  NJH'iER  2 

BEGINNING  WAVENUMBER  OF  645.300  ANO  ENDING  WAVENUMBER  OF  692.700 


VEL 

PRESSURE 

TEMPEPATURE 

AVERAGED  TRANSMITTANCE 

sssassassssx 

= 33z3;:xs  = = = = = ’-zs=s: 

SS3S3SSSX3S3Sr-.  ssx  SSS'SS  T*: 

B S S S S S' 

33 

0.3000E 

-03 

210  .0 

0. lOOOOOOOOOD 

01 

32 

0.6710E- 

-01 

218.  0 

0.10)00333300 

01 

31 

0.9510E 

00 

276.0 

0.99999999050 

00 

30 

0. 1 760E 

01 

2 70.0 

0 .99999999330 

00 

29 

0.3333E 

01 

25B.  0 

0.99999999370 

0 3 

28 

0.6520E 

01 

245.0 

0.99999997520 

00 

27 

0.  1320E 

02 

234.0 

0.99999989440 

00 

26 

0.2771E 

02 

224.0 

0. 99990951950 

00 

25 

0.3220E 

02 

223.0 

0.99999934730 

00 

24 

0. 376C5 

02 

222.0 

3.99999910580 

03 

23 

0.4370E 

02 

220.0 

0.99999878640 

00 

22 

0. 5100E 

02 

219  0 

0 .99909833920 

00 

21 

O.5950E 

02 

21.  ■ .) 

0.99999772990 

00 

20 

0.6950E 

02 

21  . ,0 

0.99999689090 

00 

19 

0. 8120E 

02 

216.0 

0.99999574100 

00 

18 

0.9500E 

02 

216.0 

0.99999415450 

00 

17 

0.1 llOS 

03 

216.0 

0.99999200650 

00 

16 

0. 13C0E 

03 

216.0 

3.99998902660 

00 

15 

0.1530E 

03 

216.0 

0.99998479750 

00 

14 

0.  1 /90E 

03 

216.0 

0.99997920050 

00 

13 

0.2090E 

03 

222.  0 

3.99997176490 

00 

12 

0.2430E 

03 

229.0 

0.99996221840 

00 

11 

0.  2810E 

03 

235.0 

0.99993019060 

00 

10 

0.3240E 

33 

242.  0 

0.99993492340 

00 

9 

0. 3720E 

03 

2 43.0 

0.99991587560 

00 

8 

0. 4260? 

03 

255.0 

3.99939234760 

30 

7 

0.4070E 

03 

261.0 

0.9993622175D 

00 

6 

0. 5540? 

03 

267.0 

0.99982604000 

03 

5 

0.623 J? 

03 

273.  0 

0.999732201 00 

30 

4 

0.7100? 

03 

279.0 

0.99972921000 

00 

3 

0.  8020? 

03 

285.0 

0.99966468720 

30  . 

2 

0.9020E 

03 

290.0 

0.99958893910 

00 

1 

0.1013E 

04 

294.0 

0.99949849900 

00 

CHANNEL  NiJMSER  3 

•BEGINNING  WAVFN'JMBFR  OF  677. AOO  ANO  ENDING  WAVENUMBER  OF  702.600 


VFL 

PRESSURE 

TEMPFRATURF 

AVERAGED  TRANSMITTANCE 

svssassx 

SSS«373SSSS3- 

= = = *=  = = 

8sa*sjsssa3  = 3s3:*=as=sas3«ass 

a s««  s«8a.sas7X3rrsssz: 

S S S S S 3 

33 

0.  33  3 3E- 

-0  3 

210.0 

3. 10333330300 

01 

32 

0.671 JE 

-01 

218.0 

0.  100  30000000 

01 

31 

0.  9510E 

00 

276,0 

0,99999999090 

00 

33 

J.1763E 

01 

2 70.0 

0. 99999099970 

00 

29 

0.3350c 

01 

253.0 

0.99999999910 

00 

28 

C.6520E 

01 

265.  0 

3.99999999670 

DO 

27 

0.1320E 

02 

236.0 

0.999999986  70 

00 

26 

0. 2770E 

02 

226.0 

0.99990996310 

00 

25 

3. 3223E 

02 

223.0 

3.99999992350 

00 

26 

0.3760F 

02 

222.0 

0.99999989620 

00 

23 

0.  6370F 

02 

220.0 

0.99999936350 

30 

22 

3.5103E 

02 

219.0 

C. 99999931090 

00 

21 

0.59506 

02 

213.0 

0.99999076370 

00 

20 

0.69536 

02 

217.  0 

0.99999965170 

00 

19 

0.3120E 

02 

216.0 

0.99999952650 

00 

18 

0.  9500E 

02 

216.0 

0.99999935610 

00 

17 

3. 1113= 

03 

216,0 

0.99999912350 

00 

16 

0. 1300E 

03 

216.0 

0.99999879610 

00 

15 

0. 1530E 

03 

216.0 

0.99999333520 

00 

16 

3.1790E 

03 

216.0 

0.99999772660 

00 

13 

0.2090E 

03 

222.0 

0.99999  639360 

00 

12 

C. 2633E 

03 

229.0 

3.99999577330 

00 

11 

0.2310E 

03 

235.0 

0.99999631100 

00 

10 

0.  32606 

03 

262.0 

0 ,99999236630 

00 

9 

3.3720E 

03 

263.0 

0.99998933320 

00 

8 

0.6260E 

03 

255.0 

0.99998653860 

00 

7 

0.  687CE 

03 

261.0 

3.99993222950 

00 

6 

0.5560E 

03 

267.0 

0.99997673960 

00 

5 

0.6280E 

03 

273.0 

0.99996992060 

00 

6 

0.71036 

03 

279.  0 

3.99996  125320 

00 

3 

0.8020= 

03 

285.0 

0.99995023610 

00  • 

2 

0.93206 

03 

290.0 

0.99993671550 

00 

1 

0.  1013E 

06 

296.0 

0.99991989680 

00 

CHANNEL  \ 

REGINNINO  WAVEN'JMBF.R  OF  686.100  ANO  ENDING  WAV^NUMBCR  OF  717,900 
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VEL 

PRESSURE 

TEMPERATURE 

AVERAGED  TRANSMITTANCE 

:S3S3S3XS 

II 

II 
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=S==3S=SSS3S3==S==aZZX 

XSSSX3S3X3X53XS3X'S  SS3S3XI 

B S S S «S« 

33 

0.3000E-03 

210  .0 

0 .10000009000 

01 

32 

0.67 lOE-0  1 

218.0 

1.19900009390 

91 

31 

0.9510E 

00 

276.0 

0. 99909999980 

00 

30 

0. 1760E 

01 

270.0 

0 .99999999940 

00 

29 

0.3330F 

01 

258.  0 

0.99999999770 

00 

29 

0.6520F 

01 

245.0 

0.99999  999100 

00 

27 

0.  1320E 

02 

234.0 

1.99999996600 

00 

26 

0.2770E 

02 

224.0 

0.99999990110 

00 

25 

0.3220E 

02 

223.0 

0 .99999986200 

00 

24 

0.3760E 

02 

222.  0 

1.99999935370 

99 

23 

0.4370E 

02 

220.0 

0.99999963210 

09 

22 

0. 510CF 

02 

219.0 

0 .99999930900 

00 

21 

0.5950E 

02 

218.  0 

0.99999976200 

00 

20 

0.69506 

02 

217.0 

0.99999571610 

00 

19 

0.  3120E 

02 

216.0 

0.99999966350 

10 

18 

0.9500E 

02 

216.0 

0. 999  99959140 

00 

17 

0.1110= 

03 

216.0 

0.99999950540 

00 

16 

0. 1399E 

0 3 

216.  0 

9.59999939570 

00 

15 

0 .15306 

03 

216.0 

0.99999925310 

00 

14 

0.  1750E 

03 

216.0 

0 .59999916  160 

00 

13 

0. 2090E 

03 

222.0 

0.999  )9885010 

00 

12 

0.2430E 

03 

229  .0 

0.99999352350 

00 

11 

0.  2310E 

03 

235.0 

9.99999317490 

90 

10 

0.3240E 

03 

242.0 

0.99999745310 

00 

9 

0.3720E 

03 

2 48.0 

0.99999663170 

00 

8 

0.4260E 

03 

255.  0 

1.99999552470 

09 

7 

0.4870E 

03 

261.0 

0.99999404810 

00 

6 

0.  5540E 

03 

267.0 

9 .99999216950 

00 

5 

0.6230E 

03 

273.0 

0.99993974680 

00 

4 

0.7100E 

03 

279.0 

0.99993666430 

00 

3 

0.  8020E 

03 

285.0 

0.99998269350 

90 

2 

0.90206 

03 

290.0 

0.99997773920 

00 

1 

0.10136 

04 

294.0 

0.99997168330 

00 

CHANNEL  NUMBER  5 

beginning;  wavenumber  of  696. SOO  and  ending  wavenumber  of  733.500 


LEVcl 

PRESSURE 

TEMPERATURE 

AVERAGED  TPANSMI TTA 

33 

3.  30.33E 

-03' 

210.0 

0.  ICO  OOOOOOOD 

01 

32 

0.6710E' 

-01 

218.0 

0. 10003000000 

01 

31 

0.  9510n 

00 

276.0 

0.1'30030J  3 330 

31 

30 

J.1760E 

01 

270.0 

0. lOOOOCOOOOO 

01 

29 

0.3330E 

01 

258.0 

0.99999999990 

00 

28 

0.6520E 

01 

245.0 

0.99999999940 

03 

27 

0.1320E 

02 

234.0 

0.99999999740 

00 

26 

0. 277CE 

02 

224.0 

0 .99999999120 

00 

25 

0. 32200 

0 2 

223.0 

0.99999998920 

00 

24 

0.3760E 

02 

222.0 

0.99999998660 

00 

23 

0.  43  70F 

02 

220.0 

0.99999998350 

30 

22 

J.5103E 

02 

219.0 

0.99999997970 

00 

21 

0.5950E 

02 

218.0 

0.99999997530 

00 

20 

0.6950E 

02 

217.0 

0.99999997330 

0) 

19 

0.8120E 

02 

216.0 

0.99999996490 

30 

18 

0.  9500E 

02 

216.0 

0.99999995880 

00 

17 

).  line 

03 

216.0 

0.99999995210 

00 

16 

0.13005 

03 

216.0 

0.99999994410 

00 

15 

0.  15305 

03 

216.0 

3.99999993410 

30 

14 

0.1790E 

03 

216.0 

0.99999992220 

00 

13 

0.20905 

03 

222.0 

0 .99999 9';'0690 

00 

12 

0.24306 

03 

229.0 

3.99999983600 

00 

11 

0.28105 

03 

235.0 

0.99999985760 

00 

10 

0.  32405 

03 

242.0 

0.99999981810 

00 

9 

0.3720E 

03 

248.0 

0.99999976410 

00 

8 

0.4260E 

03 

255.0 

0.99999968980 

00 

7 

0. 48705 

03 

261.0 

0.99999958770 

30 

6 

0.5540E 

03 

267.0 

0.99999945300 

00 

5 

0.62805 

03 

273.0 

0.99999927530 

00 

4 

0.710JE 

03 

279.0 

3.99999904120 

00 

3 

0.80206 

03 

285.0 

0.99999873020 

00 

2 

0.90206 

03 

290.0 

0.99999833420 

00 

1 

0.10136 

04 

294.0 

0.99999782750 

00 

CHANNEL  NUHIEt  6 

(BEGINNING  WAVENiiMDFP  OF  719.^00  AND  ENOING  WAVENUN8E(^  OF  790.600 
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LEVEL 

09ESSU9E 

TEM^caATURE 

AVE^AGEO  TRANSMITTANC6 
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33 

0. 30005 

-03 

210.0 

0.10300900000 

01 

32 

0.471  )S' 

-01 

210.  0 

3.10131011130 

31 

31 

0.05105 

00 

276.0 

0.10033030100 

01 

30 

0. 1760= 

01 

270.0 

0.10030033000 

01 

29 

1.33305 

01 

253.0 

0.99999999990 

00 

28 

0.6520E 

01 

245.0 

0.99999999970 

00 

27 

0. 13205 

02 

234.0 

0.99991999930 

01 

26 

0.27706 

02 

224.0 

0.99999999610 

00 

25 

0.32205 

02 

223.0 

0.99999999490 

00  . 

24 

1.3761E 

02 

222.0 

0,99999999320 

00 

23 

0.43705 

02 

220.0 

0.99999999110 

00 

22 

0. 51005 

02 

219.0 

3.99999998320 

00 

21 

3.59505 

02 

218.0 

0,99999998440 

00 

20 

0.6953E 

02 

217.0 

0.99999991920 

00 

19 

0.81205 

02 

216.0 

0,99999997240 

00 

18 

0.95005 

02 

216.0 

0.99991996300 

00 

17 

0.11106 

03 

216.0 

0.99999995040 

00 

16 

0.13035 

03 

216.  0 

1.99999993290 

13 

15 

0. 15305 

03 

216.0 

0,99999990310 

00 

14 

0. 17906 

03 

216.0 

0.99999987540 

00 

13 

0.20916 

03 

222.0 

0.99999982390 

00 

12 

0.24306 

03 

229^0 

0.99999976130 

00 

11 

0.23105 

03 

235.0 

3.99999966530 

33 

10 

0.32405 

03 

242.0 

0.99999952833 

00 

9 

0.37205 

03 

248.0 

0.99999933950 

00 

8 

1.42616 

03 

255.0 

1.99999917720 

JO 

7 

0.48705 

03 

261.0 

0.99999671580 

00 

6 

0. 55406 

03 

267.0 

0.99999323380 

00 

5 

1.62816 

03 

273.0 

0.99999761220 

00 

4 

0.71006 

03 

279  .0 

0.99999679400 

00 

3 

0.80205 

03 

285.0 

0.99999572140 

30 

2 

0.90205 

03 

290.0 

0.99999437800 

00 

1 

0. 10136 

04 

294.0 

0.99999269020 

00 
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01 

32 
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01 

31 
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CO 
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11 

30 
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01 
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01 

20 
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01 
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01 

28 

0.6521E 

0 1 

245.  0 
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01 

27 

0. 1320E 

02 

234.0 

0.9999999<<'-,4D 

00 

26 

0.  ■’7  7 )E 

02 

224.0 

0.99999999830 

00 

25 

0.322)E 

02 

223.0 

0.99999999790 

00 

2A 

0.3760E 

02 

222.0 

0.99999999750 

00 

23 

0. 6370E 

02 

220.0 

1.99999999710 

01 

22 

0.510aE 

02 

219.0 

J.  99999999651 

00 

21 

0.  5950E 

02 

218.0 

0 .99999  9995  91 

00 

20 

1.  o951E 

02 

217.  0 

1.99999999520 

11 

19 

0.3120E 

02 

216.0 

0.99999999461 

00 

18 

0.  9500E 

02 

2 U> . 0 

3.99999999391 

00 

17 

0.11  HE 

03 

216.0 

3. '■•9991)  /933.1 

.10 

16 

0.1300C 

03 

216  -0 

0.<'y9  ‘.99992  5C 

00 

15 

0.  1530P 

03 

216.0 

1.9999999910-1 

10 

lA 

0. 17905 

03 

2:6.0 

0.99999999150 

00 

13 

0.  20  90E 

03 

222  .0 

1 .99999993920 

00 

12 

1.  24  3 JE 

03 

229.  0 

1.99999993720 

01 

11 

0.28105 

03 

235.0 

1.99999998430 

00 

10 

0.  32405 

03 

242  .0 

3.99999  99<110C 

30 

9 

1.37205 

03 

243.  0 

0.99999997371 

0.1 

3 

0.4260c 

03 

255.0 

0.99999996450 

00 

7 

0.  4370- 

03 

261.0 

1.999999911 10 

11 

6 

0.55405 

03 

267.0 

0.99999993230 

00 

5 

0.  62305 

03 

273.0 

0.9999999062.1 

00 

4 

0.710  06 

03 

2 79.0 

1.99999937  110 

11 

3 

0.80205 

03 

235.0 
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APPENDIX  H 

AVERAGED  TRANSMITTANCE 


The  result  of  the  averaged  transmittance  for  the  first  HIRS 
channel  la  listed.  Ten  monochromatic  transmlttances  were  averaged 
over  the  band.  Each  separate  transmission  calculation  Included  only 
lines  that  were  within  a distance  of  5 wavenumbers  of  the  transmission 
frequency . 
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CHANNEL  NllMaER  I 
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